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A  class  of  statistics  is  proposed  for  the  problem  of  testing  for 
location  difference  using  censored  matched  pair  data.  The  class 
consists  of  linear  combinations  of  two  conditionally  independent 
statistics  where  the  conditioning  is  on  the  number,  N  ,  of  pairs  in 
which  both  members  are  uncensored  and  the  number,  N„,  of  pairs  in  which 
exactly  one  member  is  uncensored.  Since  every  member  of  the  class  is 
conditionally  distribution-free  under  the  null  hypothesis,  H  :  no 
location  difference,  the  statistics  in  the  proposed  class  can  be 
utilized  to  provide  an  exact  conditional  test  of  H  for  all  N.  and  N.. 

If  n  denotes  the  total  number  of  pairs,  then  under  suitable 
conditions  the  proposed  test  statistics  are  shown  to  have  asjnnptotlc 
normal  distributions  as  n  tends  to  infinity.  As  a  result,  large  sample 
tests  can  be  performed  using  any  member  of  the  proposed  class. 


A  method  that  can  be  used  to  choose  one  test  statistic  from  the 
proposed  class  of  test  statistics  is  outlined.  However,  the  resulting 
test  statistic  depends  on  the  underlying  distributional  forms  of  the 
populations  from  which  the  bivariate  data  and  censoring  variables  are 
sampled. 

Simulation  results  indicate  that  the  powers  of  certain  members  in 
the  class  are  as  good  as  and,  in  some  cases,  better  than  the  power  of  a 
test  for  H  proposed  by  Woolson  and  Lachenbruch  in  their  paper  titled 
"Rank  Tests  for  Censored  Matched  Paris"  appearing  on  pages  597-606  of 
Biometrika  in  1980.  Also,  unlike  the  test  of  Woolson  and  Lachenbruch, 
the  critical  values  for  small  samples  can  be  tabulated  for  the  tests  in 
the  new  class.  Consequently,  members  of  the  new  class  of  tests  are 
recommended  for  testing  the  null  hypothesis. 
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CHAPTER  ONE 
INTRODUCTION 

1.1  What  is  Censoring? 
In  many  experiments  it  is  important  to  compare  two  or  more  treat- 
ments by  observing  the  times  to  "failure"  of  objects  subjected  to  the 
different  treatments.  For  example,  one  might  be  interested  in  com- 
paring two  drugs  to  see  if  there  is  a  difference  in  the  length  of  time 
the  drugs  (treatments)  are  effective.  The  time  period  in  which  a 
treatment  is  effective  is  called  the  survival  time  or  the  lifetime  of 
that  treatment.  Unfortunately,  in  many  experiments  the  survival  times 
for  all  treatments  under  study  are  not  observable  due  to  a  variety  of 
reasons.  If  for  a  particular  subject  the  treatment  is  still  effective 
at  the  time  of  the  termination  of  the  experiment,  then  the  correspond- 
ing survival  time  is  known  to  be  longer  than  the  observation  time  of 
that  subject.  In  that  case  the  survival  time  is  unknown  and  said  to  be 
right  censored  at  the  length  of  time  the  subject  was  under  observation. 
It  is  also  possible,  that  at  the  time  of  first  observation  of  a  treated 
subject,  the  treatment  is  already  ineffective  so  that  the  corresponding 
survival  time  is  not  observable,  but  is  known  to  be  less  than  the 
observation  time  of  that  subject.  In  that  case  the  survival  time  which 
is  unknown  is  said  to  be  left  censored  at  the  observation  time.  It 
should  be  noted  that  censored  data  can  arise  in  cases  where  responses 
are  not  measured  on  a  time-scale.  Miller  (1981)  gives  an  example  from 
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Leavitt  and  Olshen  (1974)  where  the  measured  response  is  the  amount  of 
money  paid  on  an  insurance  claim  on  a  certain  date.  For  this  example, 
ah  uncensored  response  would  be  the  amount  paid  on  the  total  claim 
while  a  right-censored  response  would  be  the  amount  paid  to  date  if  the 
total  claim  is  unknown. 

Censoring  can  arise  in  different  situations  determined  by  the 
experimental  restrictions  on  the  observation  time.  For  example.  Type  I 
censoring  occurs  when  the  period  of  observation  for  each  subject  is 
preset  to  be  a  certain  length  of  time  T;  the  survival  times  longer  than 
T  would  be  right-censored  while  those  survival  times  less  than  T  would 
be  observed  and  are  known  as  failure  times.  Type  II  censoring  occurs 
when  the  period  of  observation  will  terminate  as  soon  ^s  the  r 
failure  is  observed  where  r  is  a  predetermined  integer  greater  than  or 
equal  to  one.  If  n  is  the  total  number  of  subjects  in  the  experiment, 
then  the  n-r  survival  times  which  are  unobserved  are  right-censored  at 
the  observation  time  T.  Random  censorship  is  a  third  form  of  censoring 
and  in  many  cases  it  is  a  generalization  of  Type  I  censoring.  Random 
censoring  occurs  when  each  subject  has  its  own  period  of  observation 
which  need  not  be  the  same  for  all  subjects.  This  could  occur,  for 
instance,  when  the  duration  of  the  experiment  is  fixed,  but  the  sub- 
jects do  not  necessarily  enter  at  the  beginning  of  the  experiment.  In 
that  case  each  subject  could  have  its  own  observation  time.  When  the 
time  of  entry  is  random,  the  observation  time  will  also  be  random. 
Clearly,  if  a  subject's  survival  time  is  greater  than  the  subject's 
observation  time,  then  a  right-censored  survival  time  results;  other- 
wise, a  failure  time  is  observed  for  that  subject.  In  the  random 
censoring  case,  left-censored  survival  times  could  also  occur  if  at  the 
time  of  first  observation  of  the  subject  a  failure  has  already  occurred. 
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1.2  Censored  Matched  Pairs 

This  dissertation  Is  concerned  with  the  analysis  of  the  survival 
times  of  two  treatments  on  matched  pairs  of  subjects.  The  problem  of 
testing  for  the  difference  in  locations  of  two  treatments  using  cen- 
sored matched  pair  data  was  examined  by  Woolson  and  Lachenbruch  (1980) 
under  the  assumption  that  both  members  of  a  pair  have  equal  but  random 
observation  times  and  that  the  observation  times  for  different  pairs 
are  independent.  They  utilized  the  concept  of  the  generalized  rank 
vector  introduced  by  Kalbflelsch  and  Prentice  (1973)  and  also  used  by 
Prentice  (1978)  to  derive  a  family  of  rank  tests  for  the  matched  pairs 
problem.  The  tests  are  developed  by  imitating  the  derivation  of  the 
locally  most  powerful  (LMP)  rank  tests  for  the  uncensored  case 
(Lehmann,  1959).  In  spite  of  the  fact  that  the  Woolson-Lachenbruch 
family  of  rank  tests  (W-L  tests)  is  an  intuitively  reasonable  genera- 
lization of  the  LMP  rank  tests  for  the  uncensored  case,  the  complicated 
nature  of  these  tests  makes  it  difficult  to  investigate  their  theoreti- 
cal properties.  For  example,  it  Is  not  clear  that  these  tests  are  LMP 
in  the  censored  case.  Also,  the  asymptotic  normality  of  these  tests, 
as  claimed  by  Woolson  and  Lachenbruch,  is  strictly  valid  only  when  N. 
(the  number  of  pairs  in  which  both  members  are  uncensored)  and  N„  (the 
number  of  pairs  in  which  exactly  one  member  is  censored)  are  regarded 
as  nonrandom  and  tending  to  Infinity  simultaneously.  It  is  not  known 
whether  the  asymptotic  normality  holds  (unconditionally)  as  the  total 
sample  size  tends  to  infinity. 

For  small  sample  sizes  exact  critical  values  for  the  VJ-L  test 
statistics  could  be  based  on  their  (conditional  on  N  and  N„)  permu- 
tation distributions.  However,  the  deteinnination  of  these  critical 
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values  becomes  progressively  Impractical  when  either  N.  or  N-  is 
greater  than  8  and  the  other  is  greater  than  2.  An  example  of  the 
Woolson  and  Lachenbruch  test  statistic  is  presented  in  Chapter  Four. 

The  objective  of  this  dissertation  is  to  propose  a  class  of 
statistics  for  testing  the.  null  hypothesis,  H  :  no  difference  in 
treatment  effects.  Any  member  of  this  class  is  an  alternative  to  the 
W-L  statistic  for  the  censored  matched  pair  problem.  The  statistics  in 
the  proposed  class  are  not  only  computationally  simple,  but  under  H 
are  also  distribution-free  (conditional  on  N.  and  N^) .  Furthermorei 
the  critical  values  of  these  statistics  are  easily  tabulated. 

Chapter  Two  begins  with  the  general  setup  of  the  problem  and  some 
preliminary  results  which  are  needed  in  Chapter  Three .  The  class  of 
linear  combinations  of  two  conditionally  distribution-free  and  indepen- 
dent statistics  for  testing  H  is  proposed  in  Chapter  Three  and  it  is 
shown  that  the  asymptotic  normality  of  the  corresponding  test  statis- 
tics hold  unconditionally  as  the  total  sample  size  n  tends  to  infinity. 
Simulation  studies  are  used  in  Chapter  4  to  compare  the  test  proposed 
by  Woolson  and  Lachenbruch  (1980)  to  the  tests  proposed  in  this 
dissertation. 


CHAPTER  TWO 
TWO  CONDITIONAL  TESTS  FOR  CENSORED  MATCHED  PAIRS 

2.1  Introduction  and  Summary 

In  this  chapter  two  statistics  for  testing  the  null  hypothesis  H^, 
that  there  is  no  difference  in  treatment  effects,  using  censored ■ 
matched  pair  data  are  developed.  Either  of  these  statistics  , can  be 
utilized  to  test  for  a  shift  in  location,  and  together  they  provide 
conditionally  independent  tests  of  H  conditional  on  the  numbers  of 
observed  pairs  in  which  both  members  of  a  pair  are  uncensored  (N  )  and 
in  which  exactly  one  member  of  a  pair  is  censored  (N^) . 

The  two  statistics  are  proposed  in  Section  2.4  and  it  is  shown 
that  the  exact  conditional  distributions  as  well  as  the  conditional 
as3nnptotic  distributions  of  the  statistics  are, easily  specified.  Some 
notations  and  assumptions,  are  stated  in  Section  2.2,  while  Section  2.3 
contains  preliminary  results  needed  in  the  development  of  Section  2.4. 

2.2  Notation  and  Assumptions. 

This  section  begins  with  a  set  of  notation  which  will  be  used 
hereafter  throughout.  Let  {(Xj^°,Y^°,Cj,) ;  1=1,2,.  ..,n}  denote 
Independent  random  variables  distributed  as  (X  ,  Y  ,  C) .  .  In  addition, 
(X°,Y°)  and  C  are  referred  to  as  the  true  value  of  the  treatment 
responses  and  the  censoring  variable,  respectively.  The  actual  ob- 
served value  of  the  treatment  responses  is  (X,Y)  where  X  =  min  (X  ,C) 
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and  Y  =  min  (Y°,C).   In  addition  to  (X,Y)  one  also  observes  the  value 
of  the  random  variable  (6(X),6(Y))  where 

^^^^"   0   X=C 
and 

1   Y-Y° 
fiW=   0   Y=C   . 

Note  here  that  the  indicator  random  variable  6 (X)  (respectively,  6(Y)) 
indicates  whether  or  not  X  (respectively,  Y)  represents  a  true 
uncensored  response. 
Let 


N,  =  Z  6(X.)  6(Y  )  , 
^  i=l   ^ 


and 


N,,  =  E  6(X.)[1  -  6(Y.)]  , 
ZL  i^i   1        1 


N-p  =  I  [1  -  6(X,)]  6(Y.)  , 


N   s  N   +  N 

"2     2L   "2R  • 


N,  =  E  tl  -  6(X.)][1  -  6(Y  )] 
"   i=l 


=  n  -  N^  -  N2 


Clearly,  N  ,  N-,  and  N_  denote  the  number  of  pairs  (X^,Y^)  in  which 
both  members  are  uncensored,  exactly  one  member  is  censored,  and  both 
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members  are  censored,  respectively.  Similar  interpretations  for  N., 
and  N-„  may  be  given. 

The  test  procedures  under  consideration  here  are  based  on  the 
observed  differences  tD  =X.-Y,  ;  i=I,...,n}.  It  is  said  an  observed 
difference  D.  falls  into  category  CI  if  6(X^)6(Y^)=1,  into  category  C2 
if   6(X^)[1-6(Y^)]   +   [1-6(X^)]6(Y^)=1,   and   into   category   C3   if 

[1-6(X.)][1-6(Y,)]  =  1.   It  is  evident  that  N,  is  the  ,  number  of  pairs 

•  i  '     '  1     ■  ,  ^  '    .     ^ 

(X,,Y.)  in  category  Ck,  k=l,2,3.  Also,  as  noted  by  Woolson  and  Lachen- 

bruch  (1980)  a  D,  in  category  CI  equals  X^°-Y^°=  D^°;  furthermore,  a  D^ 

in  category  C2  either  satisfies  D^>  D^°  (i.e.,  D^°  is  left  censored  at 

D  )  if  6(Xpil-6(Y^)]  =  1  or  satisfies  D^<  D^°  (i.e.,  D^°  is  right 

censored  at  D)  if  [1-6(X.)] 6(Y.);  =  1.   The  order , relationship  between 

D,  and  D,°  for  those  D,  in  category  C3  is  unknovm. 
i      i  i 

Additional  notation  is  needed  for  the  Nj^  and  N2  pairs  (X^,  Y^)  in 
categories  CI  and  C2,  respectively.  Let  (X^^,  ^j^^) ,  i  =  1,  ...,  Nj  , 
denote  the  N.  pairs  in  category  CI  and  let  ^^2i*^2i^ '  ^  ~  ^'  "•*  ^2 
denote  the  N.  pairs  in  category  C2.  Moreover,  let  D^^  and  D2  be 
defined  as  follows: 


^li  =^li-^li   '  ^=  ' ^' 

and  (2.2.1) 

^2j  -^2j  -^2j    'J  -  1.  ....  N^. 
Now  a  set  of  assumptions,  which  are  used  later,  are  stated. 
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Assumptions 

(Al):   (X  °,  Y  °)  are  independent  and  Identically  distributed 

(i.i.d.)  as  a  continuous  bivariate  random  variable 

(X°.  Y°). 
(A2):  There  exists  a  real  number  0  such 

that  (X°,Y°+0)  is  an  exchangeable  pair. 
(A3):  C, ,  C-,...,  C  are  i.i.d.  as  a  continuous  random 

variable  C,  and  is  independent  of  (X  ,Y  ) . 
(AA):  C  is  independent  of  (X°,Y°). 
(A5):  Let  ^A')   indicate  probability  calculated  with  the 

value  of  0  satisfying  (A2) .   Then 

Pq{D^  e  C3}  s  PgCe)  =  Pq{X°  >  C,  Y°  >  C}  <  1. 

Note  that  in  view  of  (A2) ,  X°  and  Y°  +  0  have  identical  marginal 
distributions  and  that  (A5)  implies  PgfN^  =  n}  <  1.  Therefore,  under 
(A5),  with  positive  probability,  there  will  be  at  least  one  D^  In 
either  category  CI  or  C2  or  both.  If  the  probability  of  observing  a  D^ 
in  category  CI  is  positive,  then  Fq(x)  will  be  used  to  denote  the 
distribution  function  of  D.  .   In  other  words, 

Fg(x)  E  P^{D^^  .  X} 

(2.2.2) 

=  PqCDj  <  x  I  Dj  e  CI}. 

Finally,  because  a  test  of  H  :  0  =  0^  utilizing  D^,  i=l,2,...,n  is  a 
test  of  H  :0  =  0  utilizing  D . '  =  D.  -  0^  we  may,  without  loss  of 

O  1       1      o 

generality,  take  0^  =  0. 
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2.3  Preliminary  Results 

It  is  obvious  that  if  N2=N-=0  with  probability  pne,  then  the 
problem  under  consideration  is  the  usual  paired  sample  problem^  Since 
under  (A2) ,  the  D  are  symmetrically  distributed -about  0,  there  exist  a 
number  of  linear  signed  rank  tests  appropriate  for  testing  H  :  0  =  0 
(Randies  and  Wolfe,  1979)  of  which  the  well-known  Wilcoxon  sighed  rank 
test  is  an  important  example.  However,  if  N  >  0'  with  positive 
probability,  then  the  situation  is  different  because  a  D^  in  C2  will 
be  either  left  or  right  censored  at  D  and  not  all  the  D^  will  equal 
D  .  Lemmas  2.3.1  and  2.3.2  below  state  some  properties  of  the  D^ 
belonging  to  CI  and  C2. 
Lemma  2.3.1; 

Let  Pq{D^  e  CI}  =  P^  (0)  =  P0{X°  <  C,  Y°  <  C}. 
Suppose  0=0  and  assumptions  (Al) ,  (A2)  and  (AA)  hold.   Suppose 
further  that  P,  (0)  >  0,  and  set  i|;(X)  =  1  or  0  depending  on  whether 
X  >  0  or  X  <  0,  respectively. 

(a)  The  distribution  of  D  conditional  on  D  e  CI  is  symmetric 
about  0,  That  is,  for  every  real  number  a 

Pq{D^  <  -a  I  D^e  CI  } 

=  Pq{D^  ^  a  I  D^e  CI  }  . 

(b)  The  random  variables  |  D,  |   and  i|)(D^  )  are  conditionally  in- 
dependent. That  is,  for  every  real  number  a  SO  and  u,  =  0,  1 

Pq{|D^|   S  a,  ^|;(D.)  =  u  |  D^e  CI} 

,  =  Pq{|  D.  I  <  a  [   D.e  CI}  Po{i|;(D.)  =  u  |  D.e  CI}  . 
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Proof : 

To  prove  (a)  consider 

P  {D^  ^  -  a  I  Dj,e  Cl}  Pq{Dj,e  Cl} 

By  (A2),  X  °  and  Y  °  are  exchangeable,  therefore,  the  quantity  above 
equals 


{Y.O  -  X  °  2  -  a,  X  °  S  C,.  Y,"  S   C,} 
1      i         i     i   i     1 


=  Pq{D^  ^  a  I  D^  G  Cl}  Pq{D^  g  Cl}  , 

which  completes  the  proof  of  (a) . 

To  prove  (b)  consider  for  a  >  0, 

Pq{  |d^  I   g  a,*(D^)  =  1  I  D^e  Cl} 

=  Pq{0  <  D^  S  a  I  D^G  Cl}  . 
From  part  (a)  the  conditional  distribution  of  D^  is  syiranetric  about  0. 
Therefore,  the  quantity  on  the  right  hand  side  of  the  equation  equals 

J5[Pq{0  <  D^  g  a  I  D^G  Cl}  +  Pq{  -a  S  D^  <  0  |  D^g  Cl}] 

=  Ji[(Pq{D^  ^  a  I  D^G  Cl}  -  Pq{D^  <  0  {  D^G  Cl}) 

+  (Pq^D^  <  0  I  D^G  Cl}  -  Pq{D^  <  -  a  I  D^G  Cl})] 
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=  J2[Pq{D^  ^  a    I  D^e   Cl}  -  Po{D^=0    j  D^e   CI} 

-  Pq{D^<  -a  |D^e  Cl}] 

=  S5Pq{  -a  S  D^  S  a  |  D^  e  Cl}  -  )5PQ{Di  =  0  |  D^e  Cl}. 

Since  under  (Al)  (X°,  Y°)  has  continuous  bivariate  distribution,  then 
for  all  e  VAX°   -  Y°  =  0}  =  0,  and  since  D^e  Cl  implies  D^  =  I^i°»  then 
the  quantity  above  equals 

J5P  {  -a  ^  D^  S  a  I  D^e  Cl} 

"*^Vl°il^^  I  D^e  Cl}  .  (2.3.1) 

Now,  from  part  (a) 

Pq{D^>0  I  D^e  Cl}  =  Pq{iP(D^)  =  1  I  D^£  Cl}  =,  ^ 

and 

Pq{D^  ^  0  I  D^e  Cl}  =  ^Q^^CD^   )  =0  |  D^e  Cl}  =  h  ,. 

Consequently,  from  (2.3.1)  it  can  be  seen  that 
Pq{|D^  I  ^  a,  tKD^  )  =  1  j  D^e  CU 
=  ig  Pq{|D^|^  a  I  D^e  Cl} 

=  Pq{|D^(<  a  I  D^E  Cl}  Pq{i|'(D^)  =  1  J  D^e  Cl}  . 
Similarly, 

Pq{|D^|  ^  a.  i|;(D^)  =  0  |  D^e  Cl} 

=  PQ{|Dj<a     JD^e  Cl}  -  Po{|Dj^  a,   i|.(D^)   =  1    JD^E  Cl} 
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-  Po{|D^|  S  a  I  D^e  Cl}  PqC'I'CD^)  =  1  |  D^eCl} 
Pq{|D^|  S  a  |D^e  Cl}  [1  -  Po{'l'(D^)  =  1  |  V  Cl}] 
Pq{|D^|  S  a    I  D^e  Cl}  PqI^KD^)   =  0    |  D^eCl}    . 


This  completes  the  proof  of  (b) . 

Lemma  2.3.2;  Let  D. . ,  ...»  D.„  be  as  defined  in  equation  (2.2.1). 

Then  conditional  on  N.  =  n.,  i=l,2,3,  D  .  ,  ...,  D^^^j  ,  and  N^j^  are 

mutually  independent.   That  is,  for  real  numbers  d^ ,... ,d   ,  and 

'■  1 


integer  n 


2R 


Pq{Djj   <  dj.j=l,...,Nj,  and  N^^  =  n2^  I  N^=n^.  1=1,2,3} 

=  l\   {D    ^  d  I  N^  =  n^,  1=1.2,3}  Pq  (N^^^  =  n^^   |N^=n^,  1=1.2.3}. 
j=l      -^ 


Proof: 


Let 


P,(d.e)  =  P-{D  =  X°  -  Y°  S  d,  max  (X°,  Y°)  S  C}, 
1       o 

P2R(e)  =  Pq  {y°  s  c  <  x°  },  P2L(e)  =  Pq  tx°  g  c  <  y°}. 


and 


P2(e)   =  P2R(e)  +  P2j^(e)  =  PglD^  e  C2}   . 
Recall  that  from  assumption  (A5) 

P^Ce)   =  Pq  {min(X°,Y°)  >  C}  =  Pq  {D^  e  03}. 
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Also,  as  a  result  of  the  definition  of  P.(0)  in  Lemma  2.3.1 
PjCe)   =  Pj  (»,  e)  =  Pq{D^  e  CI}. 

For  arbitrary  Integers  "i*  i^o  ~  "2R  *   "^ZL  ^^^   ^3  ®"*^^  ^^^^ 
n  =  n.  +  n-  +  n-  and  for  an  arbitrary  permutation 

^'l W    S*l  '  *'"  S""2r'  S^"2R^1  *  ••" 


^n^  +  n2  •  ^nj  +  n2  +  1  '  •••*  V 


of 

(1,2,  ...  ,  n) 
under  assumptions  (Al),  (A3),  and  (A4)  the  following  Identity  results. 

0  11^    J  i  ij     1^ 

j  =  n^  +  1,  ....  n^  +  n^^   ;  6(X^  )[1  -  <5(Y^  )]  =  1, 

j  =  n^  +  n^j^  +  1,  ....  n^  +  n^; 

[1  -  6(X^  )][1  -  6(Y^  )]  =  1,  j  =  n^  +  n^  +  1,  ....  n} 


■{ 


n^Pj(d..e)L  [P2R<e>5   [^21^^^^  ^  [P3(e)]   .   (2.3.2) 

Therefore,  summing  (2.3.2)  over  the 
l"r"2R'"2L'"3/ 


-14- 
distinct  ways  of  dividing  n  subscripts  into  four  groups  of 
n, ,  n,^,  n-T ,  and  n,  subscripts  it  follows  that 

)|n^Pj(d  ,0)1  [T2^(Q)y^^[?2i^(e)]'^^[?^(Q){'' 

"l»"2R'''2L^"3/  b^K  J 

Now  since  (N.,  N  ,  N-)  has  a  trinomial  distribution  with  probabilities 
Pj,  (G),  i=l,2,3,  then 


V^i  "  "i   *   ^=1'2,3}  =1  )     [P^(0)]    Mp2(®)^^[V®^^^   • 


Therefore, 


Vn  ^  ^l"-"  ^N^^  \V«2R  =^2R  Pi  =  ^'   i=1.2.3} 


[      n  y  (n^  Pj(d  ,0) )    (P2R(e) )  ^^(^21^®^  ^  ^^(P3(0) )  ■ 

\"l'"2R'"2L'"3/  J°^ 

f     n  \  n  n-  n 

\nj,n2,n3/       (P^(0))    V(P2(0))  (PgCe))/" 


^2R<^>  r''^2L   <^>  C' 


.("2    V    [;i         VVQ>V[^2R^^> 
\n2^/      b=l         Pi(0)       J  \P2(®> 


P2(0)    J  1^2^®^       J  (2.3.3) 


which  shows  that  the  joint  conditional  distribution  of  D   ,  1^2^^  and  N2j^ 

is  a  product  of  the  marginals ,  thereby  completing  the  proof  of 

Lemma  2.3.2.  D 

Lemmas  2.3.1  and  2.3.2  will  be  utilized  in  the  next  section  to  develop 
two  test  statistics  for  testing  H  :  0  =  0. 
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2,4  Two  Tests  for  Difference  In  Location 
Lemma  2.3.1  suggests  that  one  could  use  the  observed  differences  In 
CI  to  define  a  conditional  test  (conditional  on  N  =  n.)  of  H  :  0=0. 
As  a  consequence  of  the  symmetry  of  the  conditional  null  distribution  of 
the  D..,  j=l,2,...,  n. ,  any  test  for  the  location  of  a  symmetric  popula- 
tion is  a  reasonable  candidate  for  a  test  of  H  :   0=0.  For  example, 

o 

one  might  consider  using  a  linear  signed  rank  test.   The  well-known 
Wilcoxon  signed  rank  test  (Randies  and  Wolfe,  1979,  p.  322)  is  an  obvious 
choice  due  to  its  high  efficiency  for  a  wide  variety  of  distributions 
and  the  availability  of  tables  of  critical  values  for  a  large  selection  of 
sample  sizes.  Accordingly,  consider  the  statistic  T.   =  "^ir.^^^  where 

T,  (9)  =  J:^i>   (D,,   -  0)  R*  (0) 
In     j^j    Ij        J 


with  R*  (0)  denoting  the  rank  of  j  D   -  0  |  among 

1^11-  M'    1^2  -M  l^-  ^1  • 

Clearly,  the  D.  in  C2  also  contain  information  about  H  .  Lemma 
2.3.2  implies  that  conditional  on  N  =  n.,  T!iy=   n.,  any  test  based  on  N-j, 
and  N-  will  be  Independent  of  the  test  based  on  T  .  Since,  as  will 
subsequently  be  shown,  under  H 

E{N2j^  "  ^2L  1  \  "  "i'  ^=1'2,3]  =  0, 

a  second  statistic  for  testing  H  is  T-   where 

T   =N   -N   =2N   -N 
^2n    2R   "2L   ^"2R    2  * 
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Some  properties  of  the  conditional  distributions  of  T.   and 
T„  follow  directly  from  Lemmas  2.3.1  and  2.3.2.   These  properties  are 
summarized  in  Lemma  2.4.1. 


Lemma  2.4.1;  Let  F  (.)  be  defined  as  in  equation  (2.2.2). 

Also,  let  E  (.)  and  Var  (•)  denote  the  expectation  and  variance  with 

respect  to  P  (.).  Under  (Al),  (A2) ,  (A3),  and  (A4) : 
9 

(a)  Conditional  on  N.  =  n . ,  1=1,2,3  ,  T,   and  T.  are  independent. 
11  In      2n 


^^^  ^e  f^lnhi  =  "i'  i=1.2.3]  =  y^^O)  , 


e 

and 


Var  [Tj^  I  N^  =  n^,  1=1,2,3]  =  o  (e)  +  O(nj^) 

In 


where 


iln(e)=  qi)   /  1  -  F  (-X)  dF  (x)   +n^      [l-F(O)]   , 


0^  (e)  =  4  n,    (^1)2    {  /[I  -  F^(-x)]2dF^(x) 
In        1    ^       _„      0       o 


-t  /"  1  -  Fg(-x)dF^(x)]2}   . 


^""^     \   ^^20  J  ^1  =  "l'  i=l'2,3]  =  ^^^ie) 
and 


Var^  [T   In  =  n.,  1=1,2,3]  =  o^  (e) 
0   2nl  1    1  2n 
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where 


and 


y2n(e)  =  n^  [P2(e)]"^  C^ZR^Q^'^ZL^^^ ^ 


a^  (e)  =  4  n2[P2(e)]"^[P2R(e^HP2L(9^^ 
2n 


Proof: 


The  proof  of  (a)  follows  directly  from  Lenma  2.3.2  since  T^^  is  only 
a  function  of  D.  ,  j=l,  ....  N^,  and  T2^  is  only  a  function  of  N^j^. 
The  proof  of  (b)  follows  by  noting  that 


n 

T 


J=l  -^ 

^1 
Z  Z  ^(J)       +  D     )    +  z     ^(D,.) 

1  S  1  <  J   S  n^  ^^         '^        j=l  'J 


where 


==   (2l>  "l.n^  *     "1  ^l.n^ 


-1 


1  l<i<JSnj^ 


-1   "1 


are  two  U-statistics  (Randies  and  Wolfe,  1979,  page  83).  Since  U, 

i  ,n- 

and  W,   are  unbiased  estimates  P  (D, ,  +  D, .  >  0)  and  P^(D,  >  0), 
l,n.  0  li    Ij  G  13 
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respectively,  then 


^in<«)  =  (2i>  V^ii^^i  >o>^"rV°ii  ^°^- 


Upon  noting 


^G^^i  -  ^j  >  0)  =       ^     ^0(^1  >  -  hi  \hi  =  ^>  '^  ^e^^ij  ^  ^> 


/   [1  -  Fg(-x)]dFQ(x) 

—  00 


and 


V^ij  >°>  =  1- V°>' 


the  expression  for  p.  (6)  follows. 
Now, 


<2l)'^"0"l.n^*-?  ^"0^1,n^ 


.  2  (^1)  n^  CoVgCU^^^^.  W^^^^) 


2  - 

("l  )   Var^  U,    +  n"^  Var_  W, 
^2     ^  0  l,n    .     0  l,n 
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where  p   (0)  =  Corr[U,     W,   I   is  the  correlation  coefficient 

between  U,    and  W,      . 
l,n.     .  l.Hj^ 

As  a  result  of  Lemma  A  part  (iii)  of  Serfling  (1980,  p.  183)  it 
follows  that 


-1 


Also, 


=  4n^-^  {P^(D^^  ^  D^2  >   °'  ^1  *   °13  '   °^ 


-  [Pg(Djj  +  D^2  >  0>1^>  +  0(nj"^) 


4nj^{  /  [1-  Fq(-x)]2  dF^Cx) 

—00 


-  [  /"l  -  Fq(-x)  dFgCx)]^}  +  0(nj"^) 


^"e^i,n^=  nr'^^et^^^^i)^' -  f^G  ^^^1^^  ^ 


=  0(nj"^)   . 
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Since  p   (e)  =  0(1),  then  it  follows  that 
"l 


^«e   [TiJN^  =  n^.i=1.2.3] 


-1  2  00  , 

=  A  nj        (^1)      {  /   [1  -  FqC-x)]     dF^(x) 


-   [/"[I  -  F„(-x)   dF^(x)]2}   +   ("D^   .   OdiT^) 

0  0  ^  '■ 


*  iij^  .  oCiij"') 


.2(;i)nj[Var^U,„^)''tO(n,-b)'» 


2 
ai„(e)     +     (2I  )      •  0(n^'^)   +  nj2   .  oCn^'b 


+  2(^l)n^    [0(n^"h][0(n^"^)] 


oijj(e)        +     0(n^^)   +  0(nj) 


.  0   (nJ/2) 


'ln(e)      -     0(n^^) 
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The  proof  of  (c)  follows  from  Lemma  2.3.2  when  It  is  noted  that 


V^2R  =  "2R  I  ^1  =  ^'  ^==1'2'3> 


fe'te}'"{ 


-k  "2L 
^2l(«>  \ 


P2  (e) 


which  indicates  that  the  conditional  distribution  of  H^^   is  Binomial  with 

parameters  n-  and  [P2(6)]   ^2R^®^' 

Now, 

^e  ^hn  I  ^'i  =  ^'i^i'Z.s] 


and 


=0fN2R  -  ^2L    I     ^i  "  "i'i=-l'2,3] 


ti2  [P2  (e)]"^  [P2R  <e)  -  P2L  (6)] 


=  y2n  <«> 


^"0  f^2n  I  h   =  ^'i=l'2,3] 

=  Var^  [2N2j^  -  N2  |  N^  =  n^. 1=1. 2. 3] 
=  4  n2  [P2  (e)]"^  tP2R^®^^  t^2L^®^^ 


=  alie) 


which  completes  the  proof  of  Lemma  2.4.1. 
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Note  that  if  0  =  0  then  Fq(x)  is  symmetric  by  part  (a)  of  Lemma 
2.3.1  and  as  a  result 


yj^(O)  =  (^1)   /  1  -  FQ(-x)dFQ(x)  +  nj[l  -  Fq(0)] 


=  ("D  /  Fq(x)  dFpCx)  +  (1/2)  .  n^ 


(1/2)  n,(n  -1)   /  u  du  +  (1/2)  n 
^  ^      0 


(1/4)  nj^(nj,  -1)  +  (1/2)  n^ 


n^  (n^  +  1) 


Similarly, 


o2  (0)  =  A_  (^1  )2  m/3)  _  (1/2)2} 
In     "l   2 


n. 


where 


nj(nj+l)(2nj+l) 

Q    _   

n. 


(2.4.1) 


a^       +  0(n?)  (2.4.2) 

"^1 


1  24 

is  the  null  variance  of  the  Wilcoxon  signed  rank  statistic  based  on  n^ 
observations. 
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Also,  since  under  (A2)  PoR^^^  "  ^2L^°^  "  ^^^^^  ^2^^^' 


^2n  <0^  "   "2 


^2R<°^  -  ^2l(0> 


P2(0) 


(2.4.3) 


Similarly, 


a   (0)  =  4n  .  2 

2"        ^      [P2(0)]^ 


(2.4.4) 


In  Lemma  2.4.2  it  is  shovm  that,  conditional  on  N  =  n^  and  N^  =  n2  > 

T,  and  T„  have  asymptotic  normal  distributions. 
In      2n 


Lemma  2.4.2; 

(a)  lim  PeUaj^(e)]"^(Tj^  -  ^^^(e))  S  X  |n^  =  n^, 1=1,2,3}  =  *(x) 
n.-x" 


where 


*(x)  =  (21?)'**  /  exp[-y^/2]  dy. 


(b)  lim  Pq^ [02^(6) ]'^(T2^  -  ^zn^®^^  ^  x  |n^  =  n^, 1=1,2,3}  =  *(x)  • 
n2-^ 
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Proof ; 


(a)   From  the  proof  of  Lemma  2.4.1 


^in  =  <2l  )   "l,n^  *  "1  "l.n^ 


Therefore,  it  follows  that 


"l''  <2l>"'  f^m  -  ^ln(«)l  = 


3/2.n,.-l, 


1        1         ^11 


Now, 


l,nj     1    .^1    Ij 
and  since  E  [f(D  )]^  =  PqCKj.  >  0)  is  finite,  then  by  corollary  3.2.5 
of  Randies  and  Wolfe  (1979)  W,    converges  in  quadratic  mean  to 

l.Hj, 

E«  W,    =  Po(D,,  >  0)  .   Also,  n^^^  ("l)"^   -^  0  as  n  -^. 
0  l,n.     9  Ij  J-     ^  ^ 


Therefore, 


p-lim  n^3/2  (H^j-l  tw^^^^.  E^w^^^^]  =  0 


as  n.-Hx>,  showing  that 


»l''  <2l  '"'  t'^ln  -  "ln<«l 


and 
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n,^     [U,    -  E_U-    ] 
1   .   l,n.    0  l.n.-' 


have  the  same  limiting  distribution  as  n.-no. 
By  taking 


'' ^1  ■  "1  <?"'' °>' 


in  Theorem  3.3.13  of  Randies  and  Wolfe  (1979),  it  is  seen  that 


l,n.    0  l.n^^ 


(^l)-^a,„(0) 


has  an  asymptotic  normal  distribution. 

Note  that  a\^   (9)   =  0(n^^)  so  that  Varg  [T^^^  |  N^  =  n^,  1-1,2,3]  -o^^"  (6) 


=  1  +  0(n,"b  which  implies  that. Var^  [T^^|N^=n^,  1=1,2*3]  •  o^^" 
tends  to  one  as  n.  -h». 


(6) 


Therefore, 


"l''<^)"'  [Tln->in(«>r     ''ln-^ln(^> 


^V    <^  'm^  ^In  <^> 


has  a  limiting  standard  normal  distribution  which  completes 
the  proof  of  (a). 
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The  proof  of  (b)  follows  since  the  conditional  distribution 
iy-   is  binomial  with  parameters  n„  and  [P„(9)]   P„j-(e) 


Lemma  2.4.3  follows  from  the  conditional  Independence  of  T   and  T^ 
shown  in  Lemma  2.4.1. 

Lemma  2.4.3; 

^m-^ln^e)       ;^2n-^2n(e> 

mi^fn  n  )^    "^^ a   (6)    '  ^'  — o      (6)    '  ^   I V  V  ^^'^.S} 

min(nj.n2)-^  In  2n 


*(x)*(y) 


The  conditional  Independence  of  T.   and  T.  indicates  that  two 
conditionally  Independent  tests  of  H  can  be  employed  using  the  same  data 
set.   In  the  next  chapter  a  test  for  H  based  on  a  combination  of 

,'■.::':    \  O  .    . 

T,   and  T„  is  considered. 
In      2n 


CHAPTER  THREE 

A  CLASS  OF  TESTS  FOR 
TESTING  FOR  DIFFERENCE  IN  LOCATION 


3.1   Introduction  and  Summary 

In  Chapter  Two  the  general  setup  was  given  for  the  problem  of 
testing  for  the  difference  in  locations  of  two  treatments  applied  to 
matched  pairs  of  subjects  when  random  censoring  is  present.  In 
addition,  two  conditionally  independent  test  statistics  were  proposed 
for  testing  H  :  0=0.  The  main  focus  of  this  chapter  is  on  a  class  of 
test  statistics  consisting  of  certain  linear  combinations  of  these  two 
test  statistics. 

In  Section  3.2  a  class  of  test  statistics  for  testing  H^:  6=0  is 
proposed  and  some  properties  of  the  conditional  distributions 
(conditional  on  N  and  N  )  of  the  members  in  this  class  are  stated.  In 
particular  it  is  shown  that  the  statistics  in  the  proposed  class  are 
conditionally  distribution-free.  The  asymptotic  null  distributions  of 
these  statistics  are  shown  to  be  normal  in  Section  3.3.  Section  3. A 
contains  a  worked  example  to  demonstrate  the  use  of  one  member  of  this 
class  for  testing  H  when  the  sample  size  is  small.  In  Section  3.5  a 
method  to  choose  one  test  statistic  from  the  class  of  test  statistics 
is  suggested. 
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3.2.  An  Exact  Test  for  Difference  in  Location 

In  this  section  a  class  of  exact  conditional  tests  of  H  based  on 

o 

T.  and  T-  is  proposed. 


Let 


Tj^(Np  =  [N^(Nj  +  1)(2N^  +  l)/24]"'^  [T^^  -  N^(Nj  +  l)/4] 


*        -J« 
^2n''  2''    2    2n 


and  {L  }  be  a  sequence  of  random  variables  satisfying: 


Condition  I 


(a)  L  is  only  a  function  of  N,  and  N-, 
n  i.  ^ 


(b)  0  ^  L  SI  with  probability  one, 


(c)  There  exists  a  constant  L  such  that  p-lim  L  =  L  as  n-*". 


The  proposed  test  statistic  is 


T^(N^,N2)  =  (1-L^^)^  T^*(N^)  .  L^^  T^^^)      . 


* 
Note  that  under  H  and  conditional  on  N.=n,,  i=l,2,  T.  (N  )  has  a 
o  XI         In  1 

distribution  which  is  the  same  as  the  distribution  of  a  standardized 


Wilcoxon  signed  rank  statistic  based  on  n  observations  (Randies  and 


Wolfe,  1979),  while  T.  (N.)  is  distributed  as  an  independent 
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standardlzed  Binomial  random  variable  with  parameters  n.  and  J5. 

Therefore,  T  (N, ,N_)  is  conditionally  distribution-free  under  H  and  a 
n   1   Z  o 


test  of  H  can  be  based  on  T  (N,,N-), 
o  n  i  / 


Hereafter,  for  the  sake  of  simplicity  T  (n.,n2)  will  denote  the 


statistic  T  (N, ,N-)  conditioned  on  N,=n. ,  i=l,2  .  Useful 
n  i  z  11 


properties  of  T  (n. ,n_)  are  stated  in  Lemma  3.2.1. 


Lemma  3.2.1:   Under  H  :  9=0  and  (Al) ,  (A2) ,  (A3),  and  (A5) , 

— — — — —         o 

(a)  Eq  [T^(n^,n2)]  =  0, 

(b)  Varp  [T^(nj,n2)]  =  1, 

(c)  T  (n,,n-)  is  symmetrically  distributed  about  0. 

n  1  z 


Proof: 


.  and  conditional  on  N^=n.,i=l,2,  T, 
o 

distribution  as  a  standardized  Wilcoxon  signed  rank  statistic. 


Under  H  ,  and  conditional  on  N  =n.,i=l,2,  T,   (N  )  has  the  same 
o  i  i        In   1 


and  T„   (N_)  has  the  same  distribution  as  a  standardized  Binomial 
2n   L 

random  variable.  Also,  conditional  on  N  =  n  ,  i=l,2,  L  is  nonrandom. 
Therefore, 

Eq[T*^(NP  |n^  =  n^,i=l,2,3]  =  0  , 
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Eq[T*^(N2)  I  N^  =  n^,i=l,2,3]  =  0  , 
Varp  [T*^(N^)  |n^  =  n^, 1=1,2, 3]  =  1   , 


and 


^"0^^2n^^2^  l^i  °  i^i'i=l'2,3]  =  1   . 

As  a  result  (a)  follows  directly,  and  by  Lemma  2.4.1(a)  it  is  clear 
that 

VarQ  [T^(ni.n2)]  =  (1-L^)  Var^  tT^*(Np  |  N^  =  n^,i=l,2,3] 

+  L^VarQ[T2*(N2)  |  N^=nj,  ,1=1,2,3] 

=  (1  -  L^)  .  L^, 

from  which  (b)  follows. 

(c)  It  is  known  (Randies  and  Wolfe,  1979)  that  the  Wilcoxon  signed 
rank  statistic  and  the  binomial  random  variable  with  p^   are 
symmetrically  distributed  about  their  respective  means.  Therefore,  it 
follows  that  the  conditional  null  distributions  of 

T,*(N,)  and  T„  (N.)  are  symmetric  about  0  under  H  .  Since  T  (N.)  and 
In  1       iin  2.  o         111  1 

T^  (N.)  are  conditionally  independent,  then  the  proof  of  Lemma  3.2.1  is 
2n  2 

complete.  '-' 

The  conditional  distribution  of  T  (N, ,N-)  is  a  convolution  of  two 

n  1  / 

well-known  discrete  distributions  (i.e.,  the  Wilcoxon  signed  rank  and 
Binomial  distributions).  The  critical  values  for  tests  based  on  T^(nj^.n2) 
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are  easily  obtained  for  small  and  moderate  values  of  n  ,i=l,2.  Note 
that  the  tables  given  in  the  appendix  are  applicable  to  the  test 
statistic 

for  n  =  1,...,10  and  n2  =  1,2,. ..,15  at  the  .01, .025, .05,  and  .10 

levels  of  significance.  Tables  of  critical  values  for  other  choices 

of  T  (n, ,n_)  can  be  similarly  produced  but  they  depend  on  the 
n  1  z 

choice  of  L  . 
n 

For  large  values  of  n  ,i=l,2,  a  test  of  H  can  be  based  on  the 
asymptotic  distribution  of  T  (N-.N.)  as  seen  in  Section  3.3. 


3.3  The  Asymptotic  Distribution  of  the  Test  Statistic 

In  this  section  it  will  be  shown  that  the  null  distribution  of 

T  (N, ,N„)  converges  in  distribution  to  a  standard  normal  distribution 
n  1  2 

as  n  tends  to  infinity,  thereby  showing  that  Tj^(Nj^,N2)  is 
asymptotically  distribution-free.  This  fact  can  be  utilized  to  derive 
a  large  sample  test  of  H  . 

The  main  result  of  this  section  is  Theorem  3.3.3  which  will  be 
proved  after  establishing  several  preliminary  results.  The  result 
stated  in  Lemma  3.3.1  is  a  generalization  of  Theorem  1  of  Anscombe 
(1952). 


Lemma  3.3.1:  Let  {T     }  for  n,=l,2,...,  n-=l,2,...  be  an  array  of 
n.,n_       1  I. 

random  variables  satisfying  conditions  (i)  and  (ii). 
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Condltion  (i):  There  exists  a  real  number  y»  an  array  of  positive 

numbers  {u   ,    }.  and  a  distribution  function  F(«)  such  that 
"l  "2 


lim       P{T      -  Y  ^  X  0)^  ^  }  =  F(x)  , 
mln(nj,n2)^   "l*"2         ''r"2 


at  every  continuity  point  x  of  F(«). 

Condition  (li) :   Given  any  e  >  0  and  n  >  0  there  exist 

V  =  v(e,n)  and  c  =  c(e,n)  such  that  whenever  mln  (nj,n2)  >  v  then 

P{|T  .   ,  -  T     I  <  e  to      for  all  n!  ,n'  such  that 
^1  n^,n^    "l'"2      "l'"2 


n!-n.j  <  en.,    jn'-n-l  <  cn^}  >  l-ri. 


Let  {n  }  be  an  increasing  sequence  of  positive  integers  tending  to 
infinity  and  let  {N.  }  and  {N.  }  be  random  variables  taking  on  positive 
Integer  values  such  that 

p-llm  (X.n  )~^N,  =  1 
*^      1  r    ir 

as  T-*^   for  0  <  X.  <  1,  1=1,2. 

Then  at  every  continuity  point  x  of  F(»), 

lim  P{T.,   „  -  Y  <  X  Uf,  „  T  r,  „  n}  =  F(x) 
r^    ^lr'^2r        [Xin^].[^2^^ 

where  [x]  denotes  the  greatest  integer  less  than  or  equal  to  x. 
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Proof; 


For  given  e>0  and  ri>0  let  v*  arid  c  satisfy  condition  (11).   Also, 


**  /        \  ** 

let  V   be  chosen  such  that  for  min  (Ajn  .X^n  )  >  v 


P{  JN^r  "  ^i"r  M  "   ^l"r^^^  ~  t  '    i=l'2. 


** 
Then  for  min  (X,n  tX^^   )  >  iiiax(v*»v  )  =  v 

^^   l^lr  ~  ^l"r  I  <cAjn^  arid   |N2^  -^2"rl^  cX2n^} 

=  P{   |Nj^  -  X^n^  I  <   cX^n^}   +     PfJN^^  -  X2nJ<   cX2n^} 

-     ^^I^Ir  -  Vrh  *^^i"r  °'^  l^2r  "  ^2^1^  ^^2"r> 


a-|)   +   (1  -|)   -   1 


=   1  -  Ti     .  (3.3.1) 


Let  n.   =  [X,n  ]  and  n2  =  [X2n  ]   and  define  the  events  A^,  B, 

,  D  ,  and  E  as  follows : 
■       r  r 

A:  In,  -  n,  |<cn,   ^In-  -n.  |<  cn_  andll^   ..  -  T^   „  |<  eu    „'  , 
r  I  Ir    Irl   Ir  'I  2r  2r '    2r    •  ^ir'^2r     Ir'  2r     Ir'  2r 

II  " 
<  £  u)      for  all  n. ,  1=1,2,  such  that 
„j  ,„2     "l'"2     '^r"2 


n  >  V,  and   |  n,  -  n.  I   <  en.  ,  1=1,2   , 
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°r''    Plr  -  "irl  ^  ^"1^  and  [  N^^  -  n^^  |  <  cn2^  , 


D  :  T„   „  -  y    ^  X   (»)„   „   and  A   , 
r  .  N^^,N2^        ^lr'"2r 


E  :  T-,    „  -  Y  ^  X  u    „  and  A   , 
r   N^^.  N^^         "ir'"2r 

where  A  '^   denotes  the  complement  of  A  . 

Note  that  for  given  0  such  that  min(nj^^,n2^)  >  v  ,  the  event  B  implies 

I  T  .   ,  -  T       |<  e  0)       for  all  n , ,  i=l,2,  such  that  n.  >  v 
I  n;.n'     "i^.  n2r'     "ir'^Zr         ^   ^^  ^ 

and  [n  -  n.j   <  en  ,  i=l,2.     Since  A^  is  equal  to  0^ 

and  I  T„   „   -  T       |<  e  u)        if  follows  that  B  and  0  imply  A_. 
'  Nlr'^2r    "lr'"2r'     "lr'^2r 

Hence,  for  given  r  such  that  min  (n.  ,n2^)  >,v, 

P{A  }  ^  P{B  and  Oy^ 

=  P{B}  -  P{B  and  0  *^} 

g  P{B}  -  P{0^*^}* 

By  condition  (ii) ,  P{B}  >  1-n  and  by  (3.3.1)  for  given  r  such  that 
min(n.  ,n  )  >  v,  P{0^}  >  1-n,  which  implies  P{0^*^}  <  n.  Therefore,  if 


min(nj^^,n2^)  >  v,  then 


P{A  }  >  l-2n  (3.3.2) 

r 


which  implies 

P{A^^}  <  2ti.  (3.3.3) 

Now, 

P{D-}  =  P{Tn    n   "  y  -''\      .  n   ^"^  ^r> 
'^      ■  Ir'  2r  Ir'   2r 

■^<Vv*  V--2/  V»2;--^^ 

"""ir-  V  '\r-  '■2r  '  Xr-  "2/ 

and  |t„    „   -  T     „  i  <e  ov,    „  } 


Since 


''»lr-  "2r  '  \r""2r'  /  '"-Ir-  "2r 


implies 


X  M    _   +  (T_   _   -  T„   „  )   <.(.x  +  e)< 


it  follows  that 


P(D  )  <  P{T    ^   -  Y  S  (x  +  e)a)„   _  } 
""  V'''2r  ''lr'''2r 
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Now 


,  D  and  E  are  mutually  exclusive  events  such  that  D  U  E 
*   r      r  ■'  r    r 


is  the  event  T„    „   -  y  =  x  u   ■,    ,  so  that 


P{D  U  E  }  =  P{D  }  +  P{E  } 
r    r      r        r 


P{T        -     y  ^    (x  +  e)   m  „   }  +  P{E  } 


<  P{T        -  Y  2  (x+e)  0)       }  +  P{A^^} 
"ir'^^Zr  "lr'"2r       "^ 


Therefore,  from  (3.3.3),  for  min  (^i^'^2t^   ^   ^' 


P{T„   „   -  Y  =  xu)       , 
N,  ,N„    '     n,  ,n-  } 
Ir  2r         Ir  2r 


S  P{T        -  Y  ^  (x+e)  0)     ^  }  +  2n   .  (3.3.4) 
"ir-  "2r  ^^Ir'  "2r 


Now,  for  min(n,  ,n-  )>v 
Ir  2r 


P{D  }  =  P{  T„    „   -  Y  ^  xw        and  A  } 
^Ir'  ^2r        "ir'  "2r     ^ 


''  \r'   «2r  '   %r'   ^2r  "  \r'   «2r^  "  '  ' 


xo)       +  (T      -  T„   „  )  and  A  }   . 
"lr'"2r     "lr'"2r   ^lr'^2r 


Now  since  the  event  A  implies 


"lr'"2r    "lr'"2r    ^lr'^2r 
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then  the  event 


{T       -  Y  g  xo)       -  eu   ■    and  A  } 
"lr'"2r        "lr'"2r     "ir'^Zr 


Implies  the  event 

{T        -  Y  <  xw        +  (T       -  T-,   „  )  and  A^} 
Therefore, 


P{D^}  S  P{T    „   -  Y  g  (x-e)  0)        and  A  } 


P{T        -  Y^  (x-e)  u     „   } 
"ir'  "2r  "ir'  "2r 


+  P{A^}  -  P{T    „   -  Y  S  (x-e)  0)        or  A  } 
"^  "lr'"2r  "lr'''2r     '^ 


Therefore,  from  (3.3.2),  If  min.(nj^^,n2^)  >  v,  then 


P(D  )  S  P{T       -  Y  ^  (x-e)  u   .    }  +  (l-2n)  -  1 
^  "lr*"2r  "lr'"2r 


^  P{T       -  Y  ^  (x-e)  (D   "    }  -  2n 
'^lr'"2r  "lr'"2r 


As  a  result,  for  given  r  such  that  min(n.  ,n2^)>v 


P  {T„       „       -     Y  ^  x(o  } 

"lr'"2r  "lr'"2r 


P{D^}  +  P{E^} 


S  P{D^} 
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>  P{T    „  -  Y  ^  (x-e)  a)„   ^  }  -  2n  .    (3.3.5) 


In  view  of  equations  (3.3.4)  and  (3.3.5)  it  follows  that  for  given 

r  such  that  min  (n,  ,n_  )  >  v 
Ir  2r 

n,  ,n»    '         n,  ,n_ 
Ir  2r  Ir  2r 


^  ^<^N,  .N  -  Y  ^  ^  03       } 
Ir  2r  Ir  2r 


<,     P{T    ^   -  Y  S  (x-e)  w    „  }  +  2n.       (3.3.6) 
n,  ,n_    '         n,  ,n_ 
Ir  2r  Ir  2r 


Furthermore,  condition  (i)  implies 


lim  P{T    „  -  Y  ^  (x-e)a)„   „   }  =  F(x-e) 

n,  .n„   '        n,  ,n„ 
r-*«     Ir  2r  Ir  2r 


and 


lim  P{T    „  -  Y  ^  (x+e)(o„   „   }  =  F(x+e) 

n,  .n-  n,  ,n„ 

r-H»     Ir  2r  Ir  2r 


at  all  continuity  points  (x-e)  and  (x+e)  of  F(.)  Therefore,  if  x  is  a 

continuity  point  of  F(.),  then  (x-e)  and  (x+e)  can  be  chosen  to  be 

continuity  points  of  F(.)  since  the  set  of  continuity  points  is  dense. 
Consequently, 

^^™  ^^\    ,N,  -y^     ^"n,  .n,  >  =  ^^^^   ' 
r-Hx>     Ir'  2r  Ir'  2r 

which  completes  the  proof  of  Lemma  3.3.1.  □ 

Lemma  3.3.2  below  follows  from  Lemmas  2.4.1,  2.4.3,  and  3.2.1. 
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Lenma  3.3.2:  Let  L  =  L  (N, ,N-)  be  defined  as  in  Condition  I, 
n    n  1  2 

Under  H  :   6=0,  and  (Al) ,  (A2) ,  (A3),  and  (A4)   , 


P-{T  (N,,N-)  ^  X   I  N.  =  n.,  i=l,2,3}  tends  to  *(x) 
U  n  i  z       '1    1 

as  min  (n.,n_)  •*^   subject  to  the  condition 


lim        L  (n,  ,n,)  =  L   . 
mm  (n-.n^)  •^" 


Proof; 

Lemma  2.4.3  shows  that  "I'^C^)  [^^^-^^^(6)]  and  ^2^^®^  ^''^2n"^2n^®^^ 
possess  conditional  distributions,  conditional  on  N  =n^,  i=l,2,3,  which 
converge  to  standard  normal  distributions  as  min(n.,n2)  tends  to 
infinity.  From  equations  (2.4.1)  and  (2.4.2)  it  is  seen  that 


M,„(0)  = 


n^(nj  +  1) 


and 


,        n,(n,  +  l)(2n  +  1) 

lim  of   (0) i— i =   1 

n  -^  ^"  24 


* 


Therefore  under  H  ,  T,   (N.)  has  a  conditional  distribution  which  is 
o   in   1 


* 


asymptotically  standard  normal.   Similarly,  under  H  :  6=0,  T2^(N2)  has  a 

conditional  distribution  which  is  asymptotically  standard  normal. 

As  a  consequence  of  Lemmas  2.4.3  and  3.2.1,  Lemma  3.3.2  holds.        □ 
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The  following  Theorem  contains  the  main  result  of  this  section. 
The  proof  makes  use  of  the  result  of  Sproule  (1974)  stated  as  Lemma  3.3.3. 

Lemma  3.3.3;   Suppose  that 

U  =  O"^    I     f(X.  ,  ...,X„  ) 

where  B  =  set  of  all  unordered  subsets  of  r  integers  chosen  without 
replacement  from  the  set  of  integers  {1,2,  ...,n}  is  a  U-statistic 
of  degree  r  with  a  symmetric  kernel  f(*).   Let  (n^}  be  an  increasing 
sequence  of  positive  integers  tending  to  ~  as  r  -^<»  and  {N^}  be  a  sequence 
of  random  variables  taking  on  positive  integer  valves  with  probability 
one.  If  E{f(X^,...X^)}^  <  », 

lim     Var(n^U)  =  r^  C,  >  0  , 
n         i 


n 


-K» 


and 


p-lim  n    N  =  1, 


then 


Proof! 


lim     P{(Ujj  -  EUjj  )    S  N^"*^  X    (r^^j)^}   =  *(x), 


It  should  be  noted  that  in  order  to  prove  Lemma  3.3.3,  only  condition 
C2  of  Anscombe  (1952)  needs  to  be  verified  since  condition  CI  is  valid 
under  the  hypothesis.   If  CI  and  C2  of  Anscombe  (1952)  are  valid  then 
Theorem  1  of  Anscombe  (1952),  and  thus  Lemma  3.3.3,  is  valid.   This  proof 
is  contained  in  the  proof  of  Theorem  6  of  Sproule  (1974).  D 
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Theorem  3.3.4:  Under  H  :  6=0,  and  (A1),(A2),  (A3),  (A4) ,  (A5) ,  and 

—————— ^—^—         o 


for  every  real  value  x. 


lim  Pq{T^(Nj,N2)  S  x}  =  *(x), 
n-><» 


Proof: 


Let  T       =  T  (n-.n,).   Then  Lemma  3.3.2  shows  that  {T     }  for 
nj^,  n^    n  1  ^  nj^,n2 

n.=l,2 n_=l,2,...,  satisfies  condition  (1)  of  Lemma  3.3.1  with  y=0 

and  0)      =1.  Note  that  under  H  , 
n,,  n-  o 


Pjl^(O)  =  p-lim  n"^N^  =  X^ 

Therefore,  from  assumption  (A5)  it  can  be  seen  that  X.  >  0  for  at  least 

one  i=l,2.  Since  P^(0)  =  0  implies  Pq{N^  >  0}  =  0  ,  then  T^(Nj^N2)  is 

equivalent  to  one  T^  (Nj  if  P . (0)  =  0  for  exactly  one  i=l,2.   In  that 
^  in  i      1 

case  by  straightforward  application  of  Theorem  1  of  Anscombe  (1952)  and 
Lemma  3.3.3  in  a  manner  similar  to  the  argument  which  follows  in  the  case 
where  X  >  0  for  i=l,2  ,  Theorem  3.3.4  follows.  Consequently,  it  will 
be  assumed  in  the  following  argument  that  X  >  0  for  i=l,2.  Now  all 
that  is  needed  is  to  show  condition  (ii)  of  Lemma  3.3.1  is  satisfied 
to  prove  Theorem  3.3.3. 

Lemma  2.4.2  shows  that  o~,  (9)  [T.  -p,  (6)]  has  a  limiting 

in      in  in 

standard  normal  distribution.  As  a  result  T^^^  (n^^)  has  a  limiting  stan- 
dard normal  distribution  under  H  :  6=0.  Now,  under  H  it  can  be  seen  that 


h>i^   -  (3np''C0,_„^-E3U,_„^)  .  o^(n^-S  ,  (3.3.7) 
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where  U,    Is  the  U-statistic  defined  in  the  proof  of  Lemma  2.4.1(b). 

l.Hj 

As  a  consequence  of  Lemma  3.3.3,  U,    satisfies  condition  C2  of 
Anscombe  (1952).   In  other  words,  for  given  e  >  0  and  n  >  0,  there  exists 
V,  and  d.  >  0  such  that  for  any  n.  >  v, 

P{(3np^  l"l  n'  "  "l  n  1^  ^  ^°^  ^^^  "l  ^"*^^  ^^^^   l"l  "  '^l  I  ^  ^l"l^ 

S  1  -  n  . 

As  a  result  of  (3.3.7)  it  can  be  seen  that  T^^  (n^  satisfies 

condition  C2  of  Anscombe  (1952).   In  other  words,  for  given 

e,  >  0  and  n  >  0  there  exists  v,  and  d^  >  0  such  that  for  any  n^  >  v^ 

P{  |Tj*  (np  -  Tj*(np  I  <  e  for  all  n^  such 

that|nj  -  nj<  d^^n^^}  S  1  -  n-  (3.3.8) 

Now  let  D„,  ...,  D_   be  as  defined  in  equation  (2.2.1),  let 
21,      in- 

D_°  denote  the  true  difference  corresponding  to  D.^  and  let 

I[D2.  <  D  °]  =  1  or  0  according  to  whether  D^^^  <  D^?  or  D^^  S  D^^      . 

Note  that  conditional  on  N_  =  n^, 
^2n  =  ^"2^  -  n^ 


=  2   J^  I  [D^,  <  D^,]   -  n. 


=  2n2  "2.  n^  "  "2 


where 
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«2 


is  the  U-statistic  estimator  of  P  [D^  "^  ^i  I^J  ^  ^^^ 

Now  by  a  similar  argument  used  above  for  Tj^^Xn^,  shows  that  V2,n 

satisfies  condition  C2  of  Anscombe  (1952).  Thus  under  H^,  T^^Cn^) 
satisfies  condition  C2.  In  other  words, for  given  e 2  >  0  a^d  n  >  0,  there 
exists  V2  and  d.  >  0  such  that  for  any  n^  >  V2 

^^  1^2n  ^"2^  "  '^2n^"2^l  ^  ^2  ^°^  ^^^  "^2   such  that 

I  n2  -  n^  \<  d^n^}     ^  1  -  n-  (3.3.9) 

Let  e  >  0  and  n  >  0  be  given  and  let  v^,V2.  <ii'^2  satisfy  (3.3.8) 
and  (3.3.9).  Let  v  =  max  (vj^,V2)  and  c  =  min  (dj^d2).  Consider 

T'n  .n  =  (^  -^^''  ^in^^l^  ^  "^  ''^2^'    '"^"^ 


'l'"2 


pH  T'    .      ,     -     T'  <  2  e  for  all  n'   n'       :  such  that 


'r"2 


n^  -  nj<  cn^,     |  n^  -  ^0^2  \^  ^^2^ 
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^  P{(1  -  L)^  |Ti*(n;)   -  T^*(np| 

+  L^  |T2*(np   -  T2*(n2^1^  ^   ^  ^°''  ^^^  "l'   ^2  ^"'^^  ^^^^ 
\n[  -  n^l    <  cn^,    |n^     -  n^\  <  01X2} 

S  P{(1  -  L)*^  l^inCni)   -  Tj*(np|    <  e     and 

n.,   n.   such  that|nj  -  nJ<  cnj^,[n2  -  n2|<  cn2} 
=  P{(1  -  L)^   |Ti*(nJ)   -  Ti*(nj)|    <  e     for  all  n^ 
such  that   jn.    -  n.  I  <  cn^} 

+     P{L^  I^Zn^^Z^   ~  ^2n^"2^l    ^  ^         ^""^  ^'^^  "2 
such  that  j  n2  -  112 1<  cn2} 

-     P{(1   -L)*^   |Ti*(nJ)   -  Tj*(nj)|    <   e 


such   that   |n.    -  n   j<  en.      and   |n2  -  n2|<  cn^} 
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S  P{(1  -  L)^  I  T^*(nJ)  -  T^*(npl  <.e   for  all  n^ 
such  that  jn.-  n.[  <  en.} 


+'  P{L^  |T2*(n2)  -  T2*(n2)|  <e   for  all  h2 


such  that  I  n-  -  1^2 1<  cn2} 
-  1   .  (3.3.10) 

Assume  0  <  L  <  1  for,  if  not,  then  by  straightforward  application  of 
Lemma  3.3.3  and  Anscombe's  Theorem  1  (1952) 

tX,N2)  =  (1  -  L)\;(Np  .  L^  T2*(N2) 

has  an  asymptotic  normal  distribution  as  n  -^o". 

.  Now  using  inequalities  (3.3.8)  and  (3.3.9)  and  applying  them 
to  (3.3.10)  with  e  =  min  {e^(l  -  L)  ^,  Zr^   ^}  ,  then 


'     1  <  2e  for  all  n, ,.  n. 
n;  ,n;     nj^.n-'  1   Z 


P{  T'  ,    ,  -  T 
^1  '"2 


such  that  jn'-  n.l<  en.,  \^2~   "2I   ^  '^'^2^ 


s  (1  -  n)  +  (1  -  n).  -  1  =  l-2n  • 


Therefore,  T  satisfies  condition  (ii)  of  Lemma  3.3.1  so  that  the 
Theorem  is  valid  for  T^(Ni,N2)  =  (1  -  L)^j*(Nj)  ^  "^2x1  ^^2^'  "^^ ^^^ 
that  the  theorem  is  true  if  L  is  replaced  by  L^  consider 
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=  [(1  -  V''  -  (1  -  L)""]  Ti*(Np  +  an""  - >"">  "^Z^^V-   <3.3.11) 

Since  T,   (N,)  and  T.  (N„)  converge  In  distribution  to  standard 
In   1       Zn  z   ,  _ 

■   ■   *  -   ii 

normal  random  variables ,  then  Tj^^CN^  and  T2^(N2)  are  0^(1) 
(Serfling,  1980,  p.  8).   Also, . since  p-llmL^  =  L  as  n  -h»,  by 
Condition  I,  then  [(1  -L^)^  -(1-L)^]  and  (L^^*  -L'^)  are  Op(l). 
Consequently  (3.3.11)  shows  that 


"^n  (^1  -  V  -  V^l'V  =°p^^?   ' 


thus  proving  Theorem  3. 3. A. 

As  a  consequence  of  the  results  obtained  in  Sections  3.2  and  3.3 

it  is  clear  that  a  distribution-free  test  of  H^:  6=0  could  be  based  on 

T  (N,  N„).  For  small  sample  sizes  exact  tests  utilizing  the 
n   1 ,   2 

conditional  distribution  of  T  (N, ,N„)  is  possible  and  for  large  sample 

n  i  / 

sizes  the  asymptotic  normality  of  T^CN^.N^)  can  be  used.   In  the  next 

section  the  use  of  T  (N,  N.)  is  illustrated  with  a  worked  example, 
ni,/ 
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3.4  A  Small-Sample  Example 

To  illustrate  the  use  of  T  (N, ,N_)  developed  in  Sections  3.2  and  3.3, 

n  i  z 

the  data  set  considered  by  Woolson  and  Lachenbruch  (1980),  which  is  a 
slightly  modified  data  set  of  Holt  and  Prentice  (1974),  will  be 
considered  in  this  section.  The  data  set  is  reproduced  below. 


TABLE  3.1 
DAYS  OF  SURVIVAL  OF  SKIN  GRAFTS  ON  BURN  PATEINTS 

Patient       1    2    3    4    5    6    7    8    9    10   11 


Survival  of 
close  match 
(exp(X.))      37   19   57+  93   16   22   20   18   63   29   60+ 


Survival  of 

Poor  Match 

(exp(Y,))      29   13   15   26   11   17   26   21   43   15   40 


Note.  57+  and  60+  denote  right-censored  observations  at  57  and  60. 
The  model  considered  by  Woolson  and  Lachenbruch  (1980)  is 

exp(X^°)  =  4,  V^^  W^ 

exp(Y^°)  =   V^^  W^  (3.4.1) 

where  ,},  >  0  is  an  unknown  parameter,  V^^  and  V^^  are  independent  and 
identically  distributed  nonnegative  random  variables  and  W^  is  an 
independent  nonnegative  random  variable  for  all  i.  Taking  the  natural 
logarithms  of  (3.4.1)  yields 
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X^°  =  0  +  log  V^^  +  log  W^ 
log  V^^  +  log  W^ 


Y.  = 

1 


(3.4.2) 


where  0  =  log  <),.  It  is  clear  that  a  test  of  H^:  ,j)  =  1  is  equivalent 
to  a  test  of  H  :  e=0.  Now  considering  X  -Y  the  data  set  is  given  in 

o  1   J- 

Table  3.2. 


TABLE  3.2 
LOGARITHMS  OF  THE  DIFFERENCE  IN  SURVIVAL  TIME  OF  THE  SKIN  GRAFTS 


Patient 


X  -  Y        0.2436   0.3795   1.3550+   1.2745   0.3747   0.2578 


Patient 


10       11 


X  -  Y        -0.2624  -0.1542  0.3819    0.6592    0.4055+ 


As  can  be  seen  from  the  data,  N^=  9,  N^j^  =  0,  ^2^7^'  ^^^  ^2^'^'     ^^ 


follows  by  some  elementary  calculations  that 


T,  =  40.  T-  =  2,  T  *(9)  =  2.073.  and  T2*(2)  =  1.414, 


In 


The  test  statistic  T  (9,2)  can  be  calculated  now  if  L  is  known, 
n  *i 

Assuming  that  nothing  more  is  known  about  the  underlying  distributions, 

a  reasonable  choice  for  L  would  be  L^  =  .5  and,  consequently, 

(1  -  L  )  =  .5.  Note  that  for  this  choice  equal  weight  is  given  to 
n 
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T,  (N,)  and  T-  (N.)  in  the  calculation  of  T  (N  N  ).   Now, 
In   1       zn   /  n   X ,  ^ 

T^(9,12)  =  (.S)*^  Tj*(9)  +  (.5)^  T2*(2) 

=  2.466. 

Comparing  this  result  with  the  critical  values  appearing  in  the 
tables  given  in  the  Appendix  it  can  be  seen  that 

P-{T  (9,2)  ^  2.298}  =  0.0093, 

On  :    J  ', 

so  that  the  observed  value  of  2.466  is  significant  at  the  .01  level. 

As  shown  in  Woolson  and  Lachenbruch  (1980)  the  W-L  test  statistic  of 

logistic  scores  yields  the  large-sample  test  statistic   Z  =  2.49, 

which  indicates  significance  at  the  .01  level.  Also,  the  exact 

level  for  the  W-L  test  statistic  based  on  the  permutation  distribution 

of  the  logistic  scores  is  11/2048  or  0.00537,  again  showing  significance 

at  the  .01  level.   Note  that  the  value  of  8/2048  as  given  in  Woolson 

and  Lachenbruch  (1980)  is  incorrect  due  to  a  minor  computational  error. 

In  the  next  section  a  method  to  choose  L  will  be  outlined. 

n   ■ 

3.5  Choosing  the  Test  Statistic 

In  section  3.4  an  example  was  presented  which  illustrated  the 

testing  of  H  :  0=0  by  using  a  linear  combination  of  Tj^(Nj^)  and 

To  (No).  In  this  section  a  method  for  selecting  the  linear  combination 
2n  2 

with  the  maximum  conditional  Pitman  Asymptotic  Relative  Efficiency 
(PARE)  is  outlined.  The  main  result  of  this  section  is  contained  in 
Theorem  3.5.3  which  Is  proved  after  establishing  some  necessary 
preliminary  results . 
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The  following  lemma  Is  a  restatement  of  the  result  appearing  in 
Callaert  and  Janssen  (1978). 
Lemma  3.5.1; 

Let  U  be  a  U-statistic  of  degree  r=2  as  defined  in  Lemma  3.3.3 
n 


such  that 


lim  Var^Cn^  U  )  =  A  C, ,  >  0. 
o     n       i 


If  U-  =  Eq   I  f(X,  X-)  j  ^  <  "   ,  then  there  exists  a  constant  M  such 
that  for  all  n^  2 


sup     I  P{n^(4Ci)"^  (U^  -  W)   ^   x}  -  ^(x)  I 


—00  <  X  ^  "' 


-3  -k 
g  M  U3  Cj  n  /^ 


Proof: 


The  proof  is  given  by  Callaert  and. Janssen  (1978) 


Lemma  3.5.2: 


Let  y     (6)   and  a^   (6)      ,   i=l,2,   be  as  defined  in  Lemma  2.4.1. 
in  in 


and 

2 


o^^B)   =Pe[Dii-^Dl2   >   0.   T>^^  +D^3>   0] 


-  ^^el°n^^i2=^°^^ 


a/(6)   =   [P2(e)]"^[P2R(e)]    [P2L(e)] 
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Suppose  that  for  some  x  >  0 


(i)  inf  0,2(6)    =  a  2     >  0   , 


(ii)  inf  0^(6)    =  0  2     >  0 

0  ^  0  ^  6     :    •  ^ 


Then 


lim     V^in'^^^^t^in^^in^^^    "^    I>1  =  "l',  ^2  "  ^2^  =  ^^'^^  ' 
n.->-" 


1=1,2  uniformly  in  0  e  [0,6]  for  every  x. 


Proof: 


The  proof  follows  from  Lemma  3.5.1  concerning  the  rate  of  convergence 
to  normality  of  the  distribution  of  a  U-statistic. 
In  Lemma  2.4.2(a)  it  was  shown  that 

.n^3'^^l)-^W,_„^-EeWj_„J   (3.5.1) 


where 


and 


l'"l    2      1  <  i  <  j  <n^     ^J     'J 


w,    =  n -^    J:^  'l'(D..) 
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It  now  will  be  shown  that  [W,     -  E  W,   ]  converges  to  0  uniformly  in 

l,nj^     0  l.n^ 

hr  ■      ,    ■■     . 

0  e  [0»:6]»  and  that  the  distribution  of  n,  [U,    -  E  U,    J  converges  to 

a  normal  distribution  uniformly  in  0  e  [0»  6]-   In  view  of  equation 

(3.5.1)  this  will  complete  the  proof  of  Lemma  3.5.2,  for  i=l.         D 

It  is  clear  that  W,    is  an  unbiased  estimator  for  P  [D,,  >  0]  and 
l.n.  0  11 

has  variance  equal  to 

-1 


Var  (W    )  =  n  ^  P^[D   >  0]  [  1  -  P  (D   >  0) ] 
9   l»nj     1    0  11  0   ii 


=  V' V°ii  >0J  V^ii^o^ 


g  (4np-^    . 
Therefore,  by  Chebyshev's  inequality  (Randies  and  Wolfe,  1979) 


Var  W, 
0  l.n, 


P  {   I  W,    -  E  W,  ^  I  >  e}  <  — r- 

0^   '  l.n^    0  l.nj^l    ^'   "  2 


2 

4  e 


which 


implies  that  [W,    -  E  W,    ]  converges  in  probability  to 
l,n^    0  l,n^ 


0  uniformly  in  0.  Since  n     (^1)~  converges  to  zero  as  n  tends 
to  infinity,  then  the  second  term  on  the  right-hand  side  of  equation 
(3.5.1)  converges  to  zero  uniformly  in  0  e[0,(S]. 
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Let  y-  =  E-[1'(D,,  +  D,.)]  .  Now,  from  Lemma  3.5.1  it  follows  that 
3    B    11     12 

for  some  positive  constant  M  and  n^  ^  2 

sup    \?^{2-'   n^*^  a^-^e)[U     -  E^U    ]  S  x  |  N^  =  n^.  N^  =  n^}  -  *(x) 

-oo<x<«>  1  1 

^  Mpg  0^(6)]"^  n^"^  .  (3.5.2) 

Note  that  u,  ^  1,  by  definition  of  the  function  "¥(•).     Further,  by 

assumption  (i)  the  right-hand  side  of  equation  (3.5.2)  is  less  than  or 

equal  to  M  a   n.  ^   for  0  S  0  g  6  .  Consequently,  the  distribution  of 

n.fU,    -  E_U,    1  converges  in  law  to  a  normal  distribution  as  n,  tends 
V    l.n^    e  l,nj'       ^  i 

to  infinity  uniformly  in  0  e  [0,6],  thus  proving  the  Lemma  for  i=l. 

The  proof  for  i=2  follows  similarly  except  that  assumption  (ii)  is 

needed  for  the  application  of  Lemma  3.5.1  to  the  U-statistic 

U„    where  U^    is  defined  as  in  the  proof  of  Lemma  3.4.1(c).       D 
2,n2       2,n2 

The  main  concern  of  this  section  is  choosing  a  test  statistic 

V  =  (1  -  ^^  ^in^^l^  "  '   <^^2>  ^'-'-'^ 

where  0  <  X  <  1.  The  next  Theorem  provides  a   method  for  selecting  X 
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such  that  the  corresponding  test  procedure  given  by 
Reject        H  if  T^^^  >  t  , 

Do  not  reiect  H     otherwise,  (3.5.4) 

•'  o 

where  t  is  a  constant,  has  the  maximum  PARE  relative  to  all  test  procedures 
corresponding  to  T, 

Theorem  3.5.3; 

2 
Let  y.  (9)  and  o^   (0),  i=l,2,  be  as  defined  in  Lemma  2.4.1. 
m        in 

Suppose  the  following  conditions  hold: 

(i)  For  some  6  >  0  and  0  e  [0,6]  , 

y.   (8)  is  differentiable  with  respect  to  0  and 
m 


p;,(o)  =  i_.  ^,(6) 


>  0,  i=l,2. 
0  =  0 


(ii)  There  exist  positive  constrants  k  ,  i=l,2, 
such  that 


k^  =  lim   [Vin^°^^"^  f^i*^  ""in^^^^'  ^"^'^• 
n  •*<*> 


(iii)   lim   tUi„(0)]'^  [''in^V^  "   ^'  ^"^'^' 
n,  •*" 

where  0  =  0(n,  ^)   . 
n      i 
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(iv)  lim     ^i„"^(0)  ^in^V  "  ^     '  ^^^'^' 
n  -^ 

where  0  =  0(n  ~^)   . 
n      1 

(v)  For  some  6  >  0,  suppose  conditions  (i)  and  (ii)  of  Lemma  3.5.1  hold. 
Then  among  all  tests  for  testing  H  :  0=0  against  H^:  0  >  0  having  the 
form  (3.5.4)  the  test  corresponding  to  the  test  statistic  T   where 

has  maximum  conditional  PARE  (as  min  (n^,  n2)  •^)  relative  to  all 

tests  of  the  form  (3.5.4). 

Proof; 

In  view  of  Theorem  10.2.3  of  Serf ling  (1980)  the  proof  is  complete 
if  the  Pitman  conditions  P(l)  through  P(6)  (see  Serf ling  (1980), 
p.  317  and  p.  323)  are  verified  for  T^^  and  T2^  conditional  on  Nj^  "  '^i' 
N,  =  n_.   Conditions  P(2)  through  P(5)  follow  directly  from  assumptions 
(i),  (ii) ,  (iii)  and  (iv),  while  P(l)  follows  from  assumption  (v)  and 
Lemma  3.5.2.  To  verify  condition  P( 6)  note  that  by  Lemma  2.4.2  the 
conditional  asymptotic  joint  distribtuion  of  (T^^,  T^^)  is  bivariate 
normal  and  from  Lemma  3.5.2  it  follows  that 

sup        sup  lPe{c,;>)[Ti, -u^,(e)]  ^x. 

OS0S6-"'<x<" 

_00  <  y  <  00 

''2n(')t^2n-^2n^'>J  =y  |Ni=nj,N2  =  n2} 

-  *(x)*(y)[   -^  0 

as  min  (n-.n-)  •*".  Thus,  the  proof  of  Theorem  3.5.3  follows  by  taking 
p  =  0  in  Theorem  10.2.3  of  Serf ling  (1980).  □ 
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Selection  of  an  optimum  test  based  on  Theorem  3.5.3  requires 
verification  of  assumptions  (1)  through  (v)  and  the  computation  of 
constants  k.  and  k..  Unfortunately,  implementation  of  this  method  of 
selection  appears  intractable  because  of  the  general  nature  of  the 
functional  form  of  F  (•)  given  in  equation  (2.2.2).   It  is  not  clear  in 
general  what  assumptions  concerning  Fp.(*)  can  be  considered  reasonable. 
However,  as  is  shown  below  in  the  special  case  where  X  and  Y  obey  the 
usual  nonparametric  analogue  of  the  paired  sample  model,  the  form  of 
F  (•)  is  considerably  simplified. 

Suppose  X  and  Y  can  be  modelled  as  follows: 


X  =  0  +  B  +  E^ 

Y°  =   B  +  E2  (3.5.5) 

where  0  is  a  constant,  B  is  a  random  variable  representing  the  "pairing" 


variable,  and  E  ,  i=l,2  are  i.i.d.  according  to  the  distribution  function 
F(»).  Furthermore,  assume  the  censoring  variable  C  has  distribution 
function  H(«)  while  C-B  has  distribution  function  G(«).  Let 


P^(x)  =  P[E^  <  (C-B)  -  X,  E^  <  (C-B)] 


=  /  F(c  -  x)  F(c)  d  G(c) 
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Note  that 


X^   =  Pq  [X°  <  C,  Y°  g  C] 

=  P  [E^  S  (C  -  B)  -  0,  E^  <    (C  -  B)] 

=  PjCe) 

and,  as  indicated  in  the  proof  of  Theorem  3.3.4,  X^,  can  be  assumed 
to  be  greater  than  0.  Therefore, 

Fg(x)  =  Pg[X°  -  Y°  S  X  I  X°  <  C,  Y°  s  C] 
A  (x,0) 


Pj(e) 


where 


A(x,  0)  =  Pq[X°  -  Y°  <  X.  X°  <  C,  Y°  <  C] 


=  P  [E^  -  E2  <  X  -  0,  E^  <  (C  -  B)  -  0,  E^  <  (C  -  B)] 

00 

=  /  P[E  -  E2  <  X  -  0,  Ej  <  c  -  0,  E2  <  c]dG(c) 
—00 

00   c 
=  /  /  P[E  ge+x-0,  E^sc-  e]dF(e)dG(c) 


00   c 
=  /  /  F[min(e  +  x,  c)  -  e]dF(e)dG(c) . 
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From  an  application  of  results  such  as  Theorems  A. 2. 3  and  A. 2. 4 

in  Randies  and  Wolfe  (1979)  it  now  follows  that  if  F(«)  is  absolutely 

continuous  and  its  derivative  f(*)  is  bounded,  then  both  Pj^(e)  and 

A(x,e)  are  dif f erentiable  functions  of  ©.   Furthermore,  Fg(«)  is  also  a 

differentiable  and  bounded  function  of  0.  Therefore,  the  representation 

for  y,  (9)  and  0^(6)  given  in  Lemma  2.4.1(b)  can  possibly  be  utilized 
In         In 

for  checking  the  desired  assumptions.   For  example,  assumptions  (iii) 
and  (iv)  follow  directly,  but  verification  of  assumptions  (i) ,  (ii) , 
and  (v)  will  require  more  conditions  on  the  specific  form  of  F(«), 
a  problem  yet  to  be  solved.   Similarly,  assumptions  (iii)  and  (iv)  are 
easily  verified  for  T-  ,  because  the  representation 

P2R(e)  =  Pe[Y°  S  C  <  X°] 

=  P[E2  S  (C  -  B)  g  E^  +  0] 

00 

=  /  [1  -  F(c  -  0)]  -  F(c)  dG(c)   , 


and 


P2L(e)  =  P0[X°  2  C  <  Y°] 


=  /  [1  -  F(c)]  -  F(c  -  0)  dG(c) 


2 
can  be  utilized  in  conjunction  with  the  formulas  for  y2n^^^  '  °2^^   given 

in  Lemma  2.4.1(c). 


CHAPTER  FOUR 
SIMULATION  RESULTS  AND  CONCLUSIONS 

A.l-  ilntrbduction  and  Sunmary 

In  the  previous  chapters  a  class  of  test  statistics  was  proposed 
for  the  purpose  of  testing  for  location  difference  in  censored  matched 
pairs.  This  class  of  test  statistics  is  the  set  of  all  standardized 
linear  combinations  of  two  conditionally  independent  test  statistics 
proposed  in  Chapter  Two.  It  was  ,  shown  in  Chapter  Three  that  under 
certain  assumptions  each  member  of  this  class  is  conditionally  distribu- 
tion-free and  has  an  asymptotically  normal  distribution. 

In  this  chapter  a  simulation  study  is  presented  in  order  to 
compare  the  powers  of  some  members  of  the  proposed  class  and  a  test 
proposed  by  Woolson  and  Lachenbruch  (1980).  In  Section  4.2  the  statis- 
tics that  are  to  be  compared  will  be  identified  and  in  Section  4.3  the 
tesults  of  the  simulation  studies  will  be  listed.  Section  4.4  will 
contain  the  conclusions  that  are  drawn  based  on  the  simulation  study. 

4.2  Test  Statistics  To  Be  Compared 
In  this  section  the  six  test  statistics  to  be  compared  using  a 
simulation  study  will  be  identified.   Four  of  these  statistics  are 
members  of  the  class  proposed  in  Chapter  Three,  another  statistic  is 
based  on  the  two  statistics  proposed  in  Chapter  Two,  while  the 
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remaining  one  is  a  member  of  the  family  of  statistics  proposed  by 
Woolson  and  Lachenbruch  (1980). 

At  this  point  some  additional  notation  will  be  needed.   Let  Z.^s 
for  i=l,2,  .  .  . ,  N.  denote  the  ordered  absolute  values  of  the 
Dj^ECi  and  let  Z^pj  =  0  and  Z^^^  +  1)  °  "*   ^^'  "i*  ^=^'^'  *  •  *  '  ^1 

denote  the  number  of  |d2.|  ,  j=l,2,  .  .  .,  ^2'   where  D2J  e  C2, 
contained  in  the  interval  [Z.^.   ^(i+1)^  while  p^,  i=  0,1,  .  .  . ,Nj 
represents  the  number  of  these  D   which  are  positive.  Also,  let 

M.  =  E  (m,  +  1)   . 
J   k=j   ^ 

The  Woolson  and  Lachenbruch  (W-L)  test  statistic  (1980)   utilizing 
logistic  scores  can  be  written  as 


T,„  =  T  +  T 
WL    u    c 


where 

N,  i 

T  =  E^[2  ^(D,.)  -  1][1  -  n  M.  /  (M  +  1)1 
u  ^^j      Ix  .^1  J     i 

N  i 

T  =  %(2p  -  m  )  +  E^2p,  -  m.)[l-J5  n  M  /(M  +1)]   . 
c      "^o    o    ^^j  "^i    1     .J  2        3 
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Under  H^,  the  variance  of  T^,  conditional  on  the  observed  pattern  of 


o 
censoring,  is  given  by 


2 


N,     i  N        i  . 

=  E^[i  -  n  M,  /  (M,  +  1)]^  +  eV[i-12  n  m  /  (m  +  i)] 

i=l     j=l  J      ^  1=1       j=l  ^      ^ 


+  1/4  m 
o 

It  should  be  noted  that  if  there  is  no  censoring,  that  is  if  N^  =  n 
with  probability  one,  then  T^  is  equivalent  to  the  Wllcoxon  signed 
rank  statistic. 

The  second  test  statistic  studied  in  Section  4.3  is  TMAX  where 

TMAX  =  MAX  {  T^*(Nj),  T2*(N2)}  • 

It  is  clear  that  TMAX  utilizes  the  test  statistics  T^^  and  T^^  which 

it  If 

are  proposed  in  Chapter  Two  since  T^^(Np  and  T2j^(Np  are  the 

standardized  versions  of  T^^  and  T2^.  Since,  according  to  Lenma 

2.4.1(a),  T   and  T.  are  conditionally  independent,  an  exact 

conditional  test  of  H  can  be  performed  using  TMAX. 
o 

The  remaining  test  statistics  are  of  the  form 

T^(N^,  Np  =  (1  -  V^  Ty^)    -  ^'^   T^^CN^)  .    (4.2.1) 
as  defined  below.  The  first  test  statistic  of  this  type  is 


T-EQ  =  (.5)^  T^*(Np  +  ).5)\*(N2) 


corresponding  to  L  =  .5  in  (4.2.1).  The  second  is 
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T-SQR  =  [N^  /  (N^  +  N^)]^  ^In^^P  ""  ^^2  '    ^^1  *   ^2^^   '^2n^^2^ 
obtained  by  selecting  L  proportion  to  N2   .  The  third  statistic  sets 

2 
L  proportional  to  N  with  the  result 


T-SS  =  [Nj  /  (N^  +  Np*S]T^^(Np  +  [N^  /  (N^  +  Np*^]  T2^(N2)   • 
The  remaining  test  statistic  to  be  studied  is  T-STD  where 


T,   +  T-  -  N,(N,+l)/4 

„  „„„         In    2n    1  1 

T-STD  =  


{[N^(Nj  +  1)(2N^  +  l)/24]  +  N2P 

Note  that  T-STD  is  obtained  by  standardizing  T^^  +  T2^  using  its  mean 

and  standard  deviation  under  H  :  6=0 

o 

Simple  calculations  can  show  that  T-STD  has  the  form  (4.2.1)  with 


K  ■  %l  '  ^%]  ^  %l^ 


where 


%^   =  Nj(N^  +  l)(2Nj+l)/24 


and 


o   =  N 
N2    2 


Therefore,  weights  given  to  T^^(Np  and  T2^(N2)  by  T-STD  are  propor- 
tional to  the  null  variances  of  T^^  and  T2^,  respectively. 
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The  statistics  T-EQ,  T-SQR,  T-SS,  and  T-STD  are  all  of  the  form 
(4.2.1)  and  are  chosen  to  represent  the  class  of  all  test  statistics 
proposed  in  Chapter  Three.    In  the  next  section  some  simulation 
results  are  given  comparing  these  statistics  among  each  other 
and  with  T-_  and  TMAX. 

WL 

4.3  Simulation  Results 
In  this  section,  following  Woolson  and  Lachenbruch  (1980),  the 
assumed  model  for  X  °  and  Y  °,  i=l,2,  .  .  n,  is  the  log  linear  model 

exp(X^°)   =  (t.  Vj^  W^ 

exp(Y^°)   =  V^^   W^ 

where  <{)  >  0  is  an  unknown  parameter,  V   and  V^^.  are  independent  and 
identically  distributed  nonnegative  random  variables,  and  W^  is  an 
Independent  nonnegative  random  variable  for  all  i.  The  simulation 
results  for  each  case  are  developed  by  generating  500  random  samples  of 
50  observations  of  (log  V   ,  log  V^j,)  and  C^  where  C^  is  an  independent 
observation  of  some  censoring  variable.   It  is  clear  that  if  e=log  ^   then, 

X^°  -  Y^°  =  0  +  (log  Vj^  -  log  V^^)    ,  i=l,2.  .  .  .,  n; 

consequently,  the  distributional  form  of  X^  -  Y^  depends  on  the 
distributional  form  of  log  V   -  log  V^^. 

Samples  from  four  distributional  forms  for  X^  -  Y^  were 
simulated  in  this  study.   The  logistic  distribution  was  simulated 
by  generating  V   and  V   with  independent  exponential  distributions. 
Another  light-tailed  distribution,  the  normal  distribution,  was 
simulated  by  generating  log  V^^  and  log  V^^  as  independent  standard 
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normal  variables .  Two  heavy-tailed  distributions ,  the  double 
exponential  distribution  and  Ramberg-Schmeiser-Tukey  (RST)  Lambda 
Distribution  (Randies  and  Wolfe,  1979,  p.  416)  were  also  studied. 
The  double  exponential  distribution  was  arrived  at  by  generating 
log  V   and  log  V   with  independent  exponential  distributions. 
The  RST  Lambda  cumulative  distribution  function  (c.d.f.)  can  not  be 
expressed  explicitly  but  its  inverse  c.d.f.  can  be  given  as  follows: 

f"^(u)  =  X^  +  [u  ^  -  (l-u)^]/A2       0  <  u  <  1. 

As  shown  in  Ramberg  and  Schmeiser  (1972)  when  X^  =0,  X3  =  ^4  =  ~  ^  ^"^ 
X  =  -3.0674,  the  RST  Lambda  Distribution  approximates  the  Cauchy 
distribution.   This  particular  choice  of  RST  Lambda  Distribution  was 
used  to. generate  log  V^.  and  log  V^^. 

To  generate  the  censoring  random  variable,  the  Uniform  [0,  B] 
distribution  was  utilized  for  the  logistic  and  double  exponential 
cases,  while  the  natural  logarithm  of  the  Uniform  [0,  B]  distribution 
was  used  in  the  normal  and  RST  Lambda  cases.  The  choice  of  B  was  made 
in  each  case  in  such  a  way  so  that,  under  H^,  the  proportion  of 
uncensored  pairs  was  approximately  75%  of  the  total  sample  size  while 
approximately  20%  of  the  total  sample  size  consisted  of  pairs  in  which 
exactly  one  member  of  a  pair  was  uncensored.  Consequently,  approxi- 
mately 5%  of  the  pairs  were  not  utilized  since  they  were  pairs  in  which 
both  members  were  censored.  Since  the  results  presented  in  this 
section  only  apply  to  the  pattern  of  censoring  described  above,  any 
conclusions  drawn  only  apply  to  this  form  of  censoring. 
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In  each  case  the  hypothesis  tested  was  H  :  0=0  against  H: .0>O,  at 
level  .05.  Based  on  the  asymptotic  distribution  of  the  test  statis- 
tics, the  critical  value  was  chosen  to  be  1.645  in  each  case  except  for 
TMAX  which  has  a  critical  value  of  1.955.  Tables  4.1,  4.2,  4.3,  and 
4.4  give  the  resulting  power  (the  proportion  of  times  the  corresponding 
test  statistic  exceeds  the  critical  value)  of  the  six  tests  described 
in  Section  4.2  for  the  logistic,  normal,  double  exponential,  and  RST 
Lambda  distributions,  respectively. 


TABLE  4.1 

POWER  OF  THE  TESTS  FOR  THE  LOGISTIC  DISTRIBUTION 

9=0     9=.181    e=.453  e=.907  .   6=1.814 

T_              .050    .186      .532  .944  1.000 
WL 

TMAX              .040    .128      .414  .900  1.000 

T-EQ              .048    .172      .512  .890  .998 

T-SQR             .054    .156      .474  .896  1.000 

T-SS              .040    .152      .444  .894  1.000 

T-STD              .040    .116      .248  .522  .922 
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T 
WL 

TMAX 

T-EQ 

T-SQR 

T-SS 

T-STD 


TABLE  4.2 
POWER  OF  THE  TESTS  FOR  THE  NORMAL  DISTRIBUTION 


6=0 

e=.141 

e=.354 

e=.707 

0=1.414 

.060 

.180 

.546 

.960 

1.000 

.046 

.142 

.344 

.886 

1.000 

.052 

.170 

.514 

.944 

1.000 

.054 

.184 

.514 

.948 

1.000 

.052 

.156 

.466 

.936 

1.000 

.050 

.150 

.418 

.876 

1.000 

TABLE  4.3 
POWER  OF  THE  TESTS  FOR  THE  DOUBLE  EXPONENTIAL  DISTRIBUTION 


9=0  e=.141  e=.354  e=.707  6=1.414 

T                .064     .250  .728  .996  1.000 
WL 

TMAX              .038     .188  .584  .970  1.000 


T-EQ 


.064 

.250 

.728 

.996 

.038 

.188 

.584 

.970 

.064 

.262 

.746 

.998 

.050 

.264 

.750 

.998 

.056 

.236 

.742 

.998 

.060 

.224 

.648 

.974 

1.000 


T-SQR  .050     .264      .750      .998       1.000 

T-SS  .056     .236      .742      .998        1.000 

T-STD  .060     .224      .648      .974       1.000 
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TABLE  4.4  . 

POWER  OF  THE  TESTS  FOR  AN  APPROXIMATE  CAUCHY  DISTRIBUTION 

6=0   e=.2    e=.4  6-=. 6  e=.8  9=1.0 

T                .050   .124     .254  .370  .528  .670 
WL 

TMAX              .048   .102     .208  .304  .470  .596 

T-EQ              .042   .126     .234  .350  .516  .652 

T-SQR             .056   .132     .288  .394  .576  .720 

T-SS              .060   .124     .272  .386  .594  .728 

T-STD              .056   .126     .248  .356  .562  .694 


Note  that  the  values  of  0  appearing  in  Tables  4.1,  4.2,  and  4.3 

2    2 
are  the  values  Oa,  .la,  .25o,  .5o,  and  lo  where  in  Table  4.1  a  =  it  /3, 

which  is  the  variance  of  the  logistic  distribution,  while  in  Tables  4.2 
and  4.3  a  =  2,  which  is  the  variance  of  the  underlying  normal  and 
double  exponential  distributions,  respectively.   In  Table  4.4  the  value 
of  e  was  chosen  such  that  the  corresponding  percentiles  of  the  RST 
Lambda  Distribution  are  approximately  equal  to  the  corresponding  per- 
centiles determined  by  the  values  of  0  used  in  Tables  4.1,  4.2,  and  4.3. 
The  simulation  results  displayed  in  this  section  show  some 
differences  in  the  powers  of  the  test  statistics  considered.   For 
example,  TMAX  appears  to  be  less  powerful  than  the  other  test 
statistics  near  the  null  hypothesis.   However,  the  results  indicate 
that  the  choice  of  critical  value  for  TMAX  is  conservative  since  the 
,05  level  of  significance  is  not  attained  under  the  null  hypothesis. 
Consequently,  if  the  critical  value  is  decreased  until  the  .05  level  of 
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slgnificance  is  attained,  then  TMAX  cam  be  expected  to  perform  better 
than  it  did  in  these  studies.  Nevertheless,  all  six  test  statistics 
are  similar  in  that  their  powers  not  only  incresse  as  0  increases,  but 
for  large  0  they  all  have  high  powers. 

In  Table  4.1  the  results  indicate  that  for  the  logistic 
distribution  T  and  T-EQ  have  similar  local  (near  the  null  hypothesis) 
powers  with  T-SQR  and  T-SS  falling  closely  behind.  Table  4.2  shows 
that  T  ,  T-EQ,  and  T-SQR  all  have  similar  local  power  when  the 
underlying  distribution  is  normal.  For  the  heavier  tailed 
distributions  the  powers  represented  in  Tables  4.3  and  4.4  show  that 
T-SQR  has  the  highest  local  power  with  the  powers  of  T^  and  T-EQ  only 
slightly  lower. 

4.4  Conclusions 
The  simulation  study  of  Section  4.3  compared  the  large-sample 
powers  for  a  test  proposed  by  Woolson  and  Lachenbruch  (1980)  and  some 
tests  from  the  class  of  tests  proposed  in  this  dissertation.  As  noted, 
in  the  previous  section  no  noticeable  difference  in  large-sample  powers 
was  observed.  Therefore,  selection  of  a  test  procedure  for  testing 
H  :  6  =  0   is  not  clear  based  on  the  results  given  in  Section  4.3. 

0 

However,  this  selection  should  not  only  take  into  account  the  power  of 
the  test  procedures,  but  also  should  be  based  on  the  availability  of 
critical  values  and  the  level  of  difficulty  involved  in  the  imple- 
mentation of  the  test  procedures. 

As  has  been  noted  in  this  dissertation,  critical  values  for  the 
proposed  test  statistics,  T  (N.,N  ),  when  the  sample  size  is  small,  are 
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easily  obtained.  This  is  not  the  case  for  the  Woolson  and  Lachenbruch 

(1980)  statistic,  T^  ,  since  the  critical  values  need  to  be  derived 

from  its  permutation  distribution  conditional  on  the  observed  scores. 

This  can  be  time  consuming  and  impractical  since  the  critical  values  of 

Tj^  depend  on  the  observed  sample  and,  thus,  can  not  be  tabulated. 

Tables  of  critical  values  for  T  (N, ,N_)  are  easily  obtained  as  noted  in 

n  1   z 

the  Appendix.  As  for  large  sample  sizes.  Theorem  3.3.4  shows  that  the 

critical  values  for  T  (N, ,N„)  can  be  found  by  using  the  standard  normal 
n  1  / 

distribution.   The  conditional  distribution  of  the  statistic  T   under 

H   is  shown  by  Woolson  and  Lachenbruch  (1980)  to  have  an  asymptotic 
o 

normal  distribution  as  min(n.,n-)  -^,  but  it  is  not  clear  if  this  holds 

unconditionally  as  the  overall  sample  size,  n,  tends  to  infinity. 

Computation  of  the  values  of  the  statistics  T^(Nj,N2)  is  not 

difficult  since  it  only  involves  the  calculation  of  the  values  of 

a  Wilcoxon  signed  rank  statistic  and  a  binomial  statistic.   The  value 

of  the  statistic  T,„  is  more  difficult  to  find  and  for  larger  sample 
WL 

sizes  is  only  easily  obtainable  through  the  use  of  a  computer. 
Therefore,  for  moderate  and  large  sample  sizes  the  lack  of  computing 
facilities  makes  the  use  of  T.^^  somewhat  impractical. 

The  simulation  study  in  Section  4.3  does  not  include  small-sample 
power  results  since  these  results  are  not  easily  obtained.  This  is  due 
to  the  difficulty  encountered  when  trying  to  find  critical  values  for 
T  ,  since  these  values  depend  on  the  particular  sample  observed. 
Consequently,  the  selection  of  the  test  procedure  can  now  be  made  based 
upon  the  availability  of  critical  values,  the  level  of  difficulty  of 
calculating  the  corresponding  test  statistic,  and  the  large-sample 
power  results  of  the  previous  section.  Based  on  these  criteria  T-EQ  is 
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the  recommended  statistic  due  to  its  ease  of  computation,  the  avail- 
ability of  tables  of  critical  values  as  in  the  Appendix,  and  since  the 
large-sample  power  results  contained  in  Section  4.3  shows  its 
performance  is  close  to  the  best  of  the  statistics  considered  in  each 
case.  However,  if  the  form  of  the  underlying  distributions  is  known, 
then  it  might  be  possible  to  apply  the  method  of  Section  3.5  to  find 

the  test  statistic  T  (N,,N_)  which  maximizes  the  conditional  Pitman 
n  1   Z 

Asymptotic  Relative  Efficiency.  In  that  case,  the  resulting  statistic 
would  be  the  recommended  choice. 


APPENDIX 


TABLES  OF  CRITICAL  VALUES  FOR 
TESTING  FOR  DIFFERENCE  IN  LOCATION 


The  tables  in  this  appendix  list  the  critical  values  of 

Tn<"l."2>  =  (•^^'^in^^)  *   (♦5)\*(n2) 

at  the  a  equal  to  .01,  .025,  .05,  and  .10  levels  of  significance  for 
n, =1,2,..., 10  and  n  =1,2,. . . ,15.  The  critical  values  for  larger  h^ 
or  n,  can  be  approximated  by  the  critical  values  obtained  from  the 
standard  normal  distribution.  When  nj=0  the  critical  values  can 
be  obtained  from  the  critical  values  of  the  standardized  Binomial 
distribution  with  parameters  n,  and  p  =  .5.  Similarly,  when  n^   =0 
the  critical  values  can  be  obtained  from  the  critical  values  of 
the  standardized  Wilcoxon  signed  rank  statistic  based  on  n^  obser- 
vations. The  critical  values  of  this  test  statistic  are  derived  for 
each  n.  and  n-  by  convolutlng  the  standardized  Wilcoxon  signed  rank 
distribution  based  on  n.  observations  with  the  standardized  Binomial 
distribution  with  n.  observations.   Since  both  distributions  are 
discrete,  exact  .01,  .025,  .05,  and  .10  level  critical  values  do  not 
always  exist,  these  levels  are  only  approximated.   In  addition  to  the 
critical  values,  the  attained  significance  level  of  each  critical  value 
(i.e.,  the  probability  that  T  (n^,n2)  is  greater  than  or  equal  to  the 
critical  value)  is  given  in  the  parentheses. 
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V 

.01 

W 

"2=2 

>2=3 

a  = 







a  = 

.025 







a  = 

.05 







a  = 

.10 





1.932      (.0625) 

.01 

n2=4 

n2=5 

n2=6 

a  = 

2.439     (.0078) 

o  = 

.025 



2.288      (.0156) 

2.439      (.0078) 

a  = 

.05 

2.121      (.0313) 

2.288      (.0156) 

1.862      (.0547) 

a  = 

.10 

2.121      (.0313) 

1.656      (.0938) 

1.862      (.0547) 

.01 

V^ 

n2=8 

n2=9 

a  = 

2.578   (.0039) 

2.707    (.0020) 

2.357    (.0098) 

a  = 

.025 

2.578   (.0039) 

2.207    (.0176) 

2.357    (.0098) 

a  = 

.05 

2.043   (.0313) 

2.207    (.0176) 

1.886    (.0449) 

a  = 

.10 

2.043   (.0313) 

1.293   (.0742) 

1.886   (.0449) 

.01 

112=10 

n2=ll 

n2=12 

a  = 

2.496   (.0054) 

2.626    (.0029) 

2.340   (.0096) 

a  = 

.025 

2.049    (.0273) 

2.200   (.0164) 

2.340   (.0096) 

a  = 

.05 

2.049   (.0273) 

1.773   (.0566) 

1.742   (.0366) 

a  = 

.10 

1.529   (.0864) 

1.638   (.0569) 

1.334   (.0985) 

.01 

n2=13 

n2=14 

112=15 

a  = 

2.472   (.0056) 

2.597    (.0032) 

2.032    (.0088) 

a  = 

.025 

1.842   (.0231) 

2.219   (.0143) 

1.985    (.0296) 

/a  = 

.05 

1.688  (.0668) 

1.561    (.0453) 

1.666   (.0299) 

a  = 

.10       , 

.1.450   (.0676) 

1.463   (.1064) 

1.301    (.0773) 

-74- 


n^=2: 

.01 

-2=1 

n2=2 

-2=3 

a  = 

-- 

— 

a  s 

.025 

.— 



• 

a  = 

.05 

■ 



2.173   (.0313) 

a  = 

.10 

1.949   (.0625) 

1.541    (.0625) 

.01 

■     h2=4 

"2=5 

n2=6 

a  = 

2.530   (.0078) 

2.681    (.0039) 

a  = 

.025 

2.363   (.0156) 

2.530   (.0078) 

2.103   (.0273) 

a  = 

.05 

1.730   (.0313) 

1.897    (.0547) 

2.048   (.0313) 

a  = 

.10 

1.656   (.0938) 

1.897    (.0547) 

1.526    (.0898) 

.01 

V^ 

n2=8 

V^ 

a  ss 

2.820   (.0020) 

2.316   (.0098) 

2.438    (.0054) 

a  = 

;025 

2.187    (.0176) 

2.316    (.0098) 

1.966    (.0273) 

a  s 

.05 

1.750   (.0586) 

/  1.684    (.0459) 

1.656    (.0688) 

a  = 

.10 

1.555    (.0742) 

1.449   (.1006) 

1.334   (.0908) 

.01 

n2=10 

n2=ll 

n2=12 

a  = 

2.290   (.0139) 

2.029   (.0098) 

2.357   (.0056) 

a  = 

.025 

1.920   (.0166) 

2.015    (.0299) 

1.949   (.0231) 

a  = 

.05 

1.658   (.0569) 

1.603   (.0380) 

1.765   (.0533) 

a  = 

.10 

1.396   (.1106) 

1.382    (.0985) 

1.919   (.1159) 

.01 

112=13 

"2=1^ 

"2=15 

a  = 

2.321    (.0120) 

2.206   (.0088) 

2.325    (.0053) 

a  = 

.025  , 

■:  2.081    (.0144) 

:       1.952   (.0243) 

1.959   (.0193) 

a  = 

.05 

1.601    (.0453) 

:     :    1.705   (.0604) 

1.594   (.0535) 

a  = 

.10 

1..297   (.1088) 

1.327    (.1229) 

1.327   (.0952) 

-75- 


.    n^=3 

.01 

;        "2=1 

"2=2 

"2=3: 

a  = 



a  = 

.025 



2.134    (.0313) 

2.359    (.0156) 

a  = 

.05 

1.756    (.0625) 

1.603    (.0469) 

a  = 

.10 

1.841    (.0625) 

1.378    (.0938) 

1.542    (.0938) 

.01 

V^ 

-2=5 

"2=6 

a  = 

2.5A8    (.0078) 

2.337    (.0078) 

2.488    (.0039) 

a  = 

.025 

2.170    (.0156) 

2.083    (.0273) 

1.911    (.0293) 

a  =. 

.05 

1.792    (.0547) 

1.705    (.0508) 

1.711    (.0625) 

a  =. 

.10 

1.414    (.1016) 

1.450    (.0977) 

1.333   (.1055) 

.01       • 

n2=7 

"2=8 

n2=9 

a  = 

2.249   (.0098) 

2.256    (.0093) 

2;312    (.0139) 

a  " 

.025     . 

2.092    (.0166) 

2.000    (.0239)    . 

1.934    (.0254) 

a  = 

.05 

1.714    (.0459) 

1.756    (.0415) 

1.650  (.0505) 

a  = 

.10 

1.401    (.1016) 

1.378    (.0908) 

1.370  (.1106) 

.01 

n2=10 

n2=ll 

n2=12 

a  = 

2.167    (.0098) 

2.248    (.0090) 

2.359    (.0120) 

a  = 

.025 

2.028    (.0299) 

1.870    (.0239) 

1.981    (.0215) 

a  = 

.05 

1.650   (.0526) 

1.773   (.0541): 

.1.603   (.0477) 

a  = 

.10 

1.272   (.1052) 

1.396    (.0917) 

1.225   (.1028) 

.01 ; 

"2=13 

n2=14 

n2=15 

d  = 

2.143   (.00.90) 

2.268    (.0124) 

2.386    (.0101) 

a  = 

.025 

1.779    (.0245) 

1.890    (.0250)     , 

2.034    (.0269) 

a  = 

.05 

1.737    (.0422); 

1:512   (.0494)  .^ 

1.682   (.0485) 

0   a 

.10 

1.408  (.0125) 

1.13A  (.0988) 

1.304   (.1128) 

-76- 


n^=4 

01 

-2=1 

"2=2 

"2=3 

a  = 

2.291    (.0156) 

2.516   (.0078) 

a  = 

025 

1.998   (.0313) 

2.033   (.0313) 

1.999   (.0234) 

a  = 

05 

1.740   (.0625) 

1.775    (.0469) 

1.741    (.0391). 

a  = 

,10 

1.481    (.0938) 

1.291    (.1094) 

1.441    (.1016) 

.01 

n2=4 

n2=5 

n2=6 

a  = 

2.189   (.0117) 

2.356    (.0059) 

2.248   (.0107) 

a  = 

.025 

1.998   (.0273) 

1.981    (.0293) 

1.929   (.0244) 

a  = 

.05 

1.740   (.0508) 

1.723   (.0430) 

1.671    (.0527) 

a  = 

.10 

1.291    (.1055) 

1.349    (.1055) 

1.352    (.0957) 

.01 

"2=7 

"2=8 

"2=5 

a  = 

2.369   (.0093) 

2.291    (.0120) 

2.379   (.0085) 

a  = 

.025 

1.871    (.0249) 

2.016    (.0251) 

1.998   (.0256) 

o  = 

.05 

1.595   (.0498) 

1.758    (.0500) 

1.695    (.0515) 

o  = 

.10 

1.318   (.1055) 

1.291    (.0994) 

1.268   (.1061) 

.01 

n2=10 

n2=ll 

1^2=12 

a  = 

2.236   (.0091) 

2.357    (.0103) 

2.299   (.0078) 

a  = 

.025 

1.927    (.0278) 

2.009    (.0211) 

1.933   (.0245) 

o  = 

.05 

.1.720   (.0475) 

1.672   (.0504) 

1.699   (.0532) 

a  = 

.10 

1.291   (.1099) 

1.312   (.0994) 

1.333  (.0991) 

.01 

112=13 

n2=14 

n2=15 

a  = 

2.281   (.0081) 

2.286   (.0101) 

2.222   (.0095) 

a  = 

.025 

1.889   (.0272) 

2.028    (.0252) 

1.946   (.0257) 

o  = 

.05 

1.631    (.0482) 

1.669   (.0518) 

1.643   (.0502) 

a  = 

.10 

1.249    (.0987) 

1.291    (.1035) 

1.322   (.1004) 

-77- 


n^=5 

.01 

"2=1 

"2=2 

"2=3 

o  = 

2.430   (.0078) 

2.274    (.0117) 

a  = 

.025 

1.947    (.0313) 

2.049   (.0234) 

1.892    (.0273) 

a  = 

.05 

1.756   (.0469) 

1.667    (.0547) 

1.701    (.0508) 

a  = 

.10 

1.375   (.1094) 

1.430   (.0938) 

1.320   (.0977) 

.01 

n2=4 

"2=5 

n.=6 

o  = 

2.272   (.0098) 

2.379   (.0098) 

2.394   (.0093) 

a  = 

.025 

2.082    (.0215) 

1.997    (.0244) 

2.008    (.0283) 

o  = 

.05 

1.700   (.0488) 

1.676   (.0498) 

1.631    (.0532) 

a  = 

.10 

1.319   (.1035) 

1.295    (.0996) 

1.255    (.0967) 

.01 

n2=7 

"2=« 

"2=9 

a  = 

2.348   (.0076) 

2.286   (.0099) 

2.317    (.0090) 

a  = 

.025 

2.004   (.0254) 

1.976    (.0220) 

1.947    (.0279) 

a  = 

.05 

1.660   (.0559) 

1.667    (.0530) 

1.666   (.0486) 

a  = 

.10 

1.316   (.0950) 

1.333   (.0988) 

1.285   (.1013) 

.01 

n2=10 

n2=ll 

n2=12 

a  = 

2.325    (.0104) 

2.351    (.0086) 

2.274    (.0099) 

a  = 

.025 

1.943   (.0252) 

1.969   (.0240) 

1.919    (.0244) 

a  = 

.05 

1.687    (.0514) 

1.643   (.0543) 

1.675   (.0520) 

a  = 

.10 

1.312   (.0940) 

1.307    (.1027) 

1.293   (.1036) 

.01 

n2=13 

ti2=14 

n2=15 

a  = 

2.242   (.0099) 

2.360   (.0095) 

2.327    (.0102) 

a  = 

.025 

2.019   (.0256) 

1.992    (.0249) 

1.962   (.0262) 

a  =» 

.05 

1.648   (.0501) 

1.617    (.0493) 

1.613   (.0526) 

a  = 

.10 

1.267    (.1028) 

1.240   (.1002) 

1.308   (.0983) 

-78- 


n.=6 

.01 

n2=l 

n2=2 

"2=3 

a  = 

2.264    (.0078) 

2.408   (.0078) 

2.337    (.0098) 

a  = 

.025 

1.967    (.0234) 

1.964   (.0273) 

1.965   (.0254) 

a  = 

.05 

1.671    (.0547) 

1.667    (.0547) 

1.668   (.0527) 

a  = 

.10 

1.374   (.1094) 

1.371    (.1016) 

1.299   (.1035) 

.01 

n2=4 

"2=5 

n2=6 

a  = 

2.264   (.0107) 

2.357    (.0098) 

2.267    (.0107) 

o  = 

.025 

1.967    (.0254) 

1.952    (.0229) 

1.970   (.0273) 

a  = 

.05 

1.671   (.0488) 

1.655   (.0498) 

1.674   (.0525) 

a  = 

.10 

1.340   (.1025) 

1.319   (.0986) 

1.260   (.1021) 

.01 

V^ 

"2=« 

"2=5 

a  = 

2.358   (.0090) 

2.315   (.0096) 

2.290   (.0102) 

a  = 

.025 

1.945   (.0236) 

1.964    (.0261) 

1.967   (.0258) 

a  = 

.05 

1.648   (.0504) 

1.667    (.0499) 

1.671    (.0518) 

a  = 

.10 

1.321    (.1010) 

1.315    (.1017) 

1.321    (.1018) 

.01 

n2=10 

n2=ll 

112=12 

a  = 

2.305   (.0099) 

2.289   (.0101) 

2.337    (.0098) 

a  = 

.025 

2.006   (.0235) 

1.975    (.0244) 

1.931    (.0245) 

a  = 

.05 

1.566   (.0503) 

1.678   (.0490) 

1.632    (.0534) 

a  = 

.10 

1.265   (.1043) 

1.307    (.1043) 

1.299   (.1023) 

.01 

n2=13 

n2=14 

n2=15 

a  = 

2.284   (.0100) 

2.313   (.0103) 

2.321    (.0096) 

a  = 

.025 

1.987    (.0244) 

1.949   (.0249) 

1.956   (.0252) 

a  = 

.05 

1.648   (.0534) 

1.653   (.0500) 

1.660  (.0505) 

a  = 

.10 

1.308   (.1000) 

1.289   (.0989) 

1.294   (.1011) 

-79- 


n^=7 

.01 

n2=l 

-2=2 

n2=3 

o  = 

2.261    (.0078) 

2.315    (.0098) 

2.300   (.0098) 

a  = 

.025 

1.902    (.0273) 

1.956    (.0273) 

1.962   (.0244) 

a  = 

.05 

1.663   (.0547) 

1.673    (.0508) 

1.703   (.0508) 

o  = 

.10 

1.305    (.1171) 

1.315    (.1055) 

1.344   (.0996) 

.01 

n2=4 

"2=5 

n_=6 

a  = 

2.261    (.0107) 

2.298   (.0095) 

2.251    (.0106) 

O    o 

.025 

2.012    (.0244) 

1.940   (.0254) 

1.971    (.0258) 

a  = 

.05 

1.663   (.0518) 

1.666    (.0513) 

1.653   (.0510) 

a  = 

.10 

1.305   (.1064) 

1.307    (.1011) 

1.294   (.1053) 

.01 

n2=7 

"2=« 

n2=9 

a  = 

2.292    (.0108) 

2.315   (.0103) 

2.261    (.0105) 

a  = 

.025 

1.941    (.0250) 

1.956    (.0265) 

1.909   (.0260) 

a  = 

.05 

1.695    (.0478) 

1.641    (.0501) 

1.663   (.0498) 

a  = 

.10 

1.287    (.0991) 

1.315    (.1025) 

1.305    (.1033) 

.01 

n2=10 

112=11 

n2=12 

a  = 

2.298   (.0103) 

2.313   (.0102) 

2.300   (.0103) 

a  = 

.025 

1.970   (.0252) 

1.954    (.0257) 

1.962   (.0249) 

a  = 

.05 

1.642   (.0522) 

1.647    (.0516) 

1.653   (.0526) 

a  = 

.10 

1.284    (.1056) 

1.305    (.1000) 

1.295    (.1034) 

.01 

n2=13 

112=14 

n2=15 

a  = 

2.295    (.0102) 

2.329    (.0096) 

2.241    (.0103) 

a  = 

.025 

1.937    (.0263) 

1.971    (.0249) 

1.982    (.0254) 

a  B 

.05 

1.664  (.0503) 

1.651    (.0499) 

1.637   (.0494) 

a  = 

.10 

1.305   (.1010) 

1.312   (.1000) 

1.278   (.0982) 

-80- 


n^=8 

.01 

"2=1 

-2=2 

"2=3 

o  = 

2.192   (.0098) 

2.287    (.0098) 

2.314   (.0093) 

a  = 

.025 

1.895    (.0273) 

1.990   (.0244) 

1.992   (.0239) 

a  = 

.05 

1.697   (.0488) 

1.693   (.0498) 

1.695   (.0488) 

a  = 

.10 

1.400   (.0947) 

1.297    (.1064) 

1.324    (.0972) 

.01 

n2=4 

"2=5 

n2=6 

o  = 

2.305   (.0098) 

2.274   (.0103) 

2.326   (.0096) 

a  = 

.025 

1.994   (.0242) 

1.977    (.0240) 

1.947    (.0269) 

a  = 

.05 

1.683   (.0505) 

1.642   (.0532) 

1.667   (.0497) 

a  = 

.10 

1.301   (.1040) 

1.306    (.1027) 

1.287    (.0997) 

.01 

n2=7 

"2=8 

"2=5 

a  = 

2.287    (.0107) 

2.287    (.0106) 

2.291    (.0100) 

a  = 

.025 

1.951    (.0253) 

1.985    (.0256) 

1.971    (.0251) 

a  = 

.05 

1.653   (.0504) 

1.688    (.0498) 

1.674   (.0498) 

a  = 

.10 

1.297    (.1032) 

1.297    (.0984) 

1.301    (.1023) 

.01 

n2=10 

n2=ll 

n2=12 

a  = 

2.284    (.0097) 

2.315    (.0097) 

2.314   (.0097) 

a  = 

.025 

1.936   (.0260) 

1.957    (.0256) 

1.930   (.0251) 

a  = 

.05 

1.639   (.0510) 

1.660   (.0500) 

1.633   (.0499) 

a  = 

.10 

1.294   (.0968) 

1.302    (.1004) 

1.299   (.1032) 

.01 

n2=13 

n2=14 

n2=15 

a  = 

2.272   (.0100) 

2.304    (.0102) 

2.299   (.0103) 

a  = 

.025 

1.971    (.0249) 

1.944    (.0259) 

1.965    (.0250) 

a  = 

.05 

1.674   (.0499) 

1.647    (.0518) 

1.643   (.0498) 

a  = 

.10 

1.285   (.0990) 

1.287   (.1007) 

1.303   (.1001) 

-81- 


^=^ 

-2=V 

n^2      ■ 

"2=3 

a  = 

01 

2.173   (.0098) 

2.298   (.0093) 

2.272   (.0110) 

a  = 

.025 

1.922   (.0244) 

1.963   (.0254) 

1.958   (.0239) 

a  = 

.05 

1.670   (.0508) 

1.634    (.0498) 

1.685   (.0496) 

a  = 

.10 

1.335   (.1064) 

1.298   (.1001) 

1.309   (.0972) 

.01 

n2=4 

n2=5 

n2=6 

a  = 

2.294   (.0103) 

2.293   (.0099) 

2.286   (.0107) 

a  = 

.025 

1.959   (.0248) 

1.958   (.0248) 

1.951    (.0262) 

o  = 

.05 

1.670   (.0505) 

1.661    (.0508) 

1.634    (.0499) 

a  = 

.10 

1;298   (.1012) 

1.288   (.1007) 

1.289    (.1038) 

.01 

"2=^ 

"2=« 

n2=9 

o  = 

2.300   (.0100) 

2.298    (.0101) 

2.309    (.0101) 

a  = 

.025 

1.965   (.0248) 

1.963    (.0251) 

1.953  (.0250). 

ot  = 

.05 

1.649   (.0508) 

1.634    (.0485) 

1.649   (.0496) 

a  = 

.10 

1.294   (.1017) 

1.296    (.1012) 

1.304   (.0997) 

.01 

n2=10 

n2=ll 

n2=12 

a  = 

2.305   (.0100) 

2.288   (.0100) 

2.282   (.0103) 

o  = 

.025 

1.943   (.0244) 

1.953   (.0242) 

1.947    (.0255) 

a  = 

.05 

1.662   (.0499) 

1.672   (.0500) 

1.634   (.0502) 

a  = 

.10 

1.300   (.0987) 

1,283   (.0980) 

1.288   (.1017) 

.01 

n2=13 

n2=14 

n2=15 

a  = 

2.283   (.0100) 

2.306    (.0100) 

2.302   (.0097) 

o  = 

.025 

1.948   (.0247) 

1.932   (.0251) 

1.960   (.0254) 

a  = 

.05 

1.639   (.0500) 

1.638   (.0500) 

1.655   (.0498) 

a  = 

.10 

1.300   (.1008) 

1.300   (.1001) 

1.290   (.1017) 

-82- 


nj=10 

.01 

n2=l 

n2=2 

n2=3 

a  = 

2.185   (.0093) 

2.261    (.0105) 

2.270   (.0107) 

a  = 

.025 

1.968   (.0210) 

1.973   (.0247) 

1.958   (.0253) 

o  = 

.05 

1.680  (.0483) 

1.685   (.0503) 

1.670   (.0505) 

a  = 

.10 

1.320   (.1079) 

1.324    (.1030) 

1.309   (.1030) 

.01 

n2=4 

n2=5 

n2=6 

a  = 

2.315   (.0097) 

2.282    (.0104) 

2.273   (.0099) 

a  = 

.025 

1.968    (.0244) 

1.938   (.0253) 

1.982   (.0249) 

a  = 

.05 

1.666   (.0508) 

1.650   (.0500) 

1.624   (.0496) 

a  = 

.10 

1.320   (.0997) 

1.289   (.1025) 

1.264   (.0997) 

.01 

n2=7 

n2=8 

n2=9 

a  = 

2.309   (.0097) 

2.324   (.0098) 

2.300   (.0103) 

a  = 

.025 

1.961   (.0244) 

1.968   (.0256) 

1.968   (.0250) 

a  = 

.05 

1.643   (.0501) 

1.676   (.0498) 

1.647    (.0507) 

a  = 

.10 

1.300   (.1005) 

1.320   (.1000) 

1.287    (.1011) 

.01 

112=10 

n2=ll 

n2=12 

a  = 

2.315   (.0097) 

2.321    (.0099) 

2.318   (.0097) 

a  = 

.025 

1.954   (.0250) 

1.967    (.0251) 

1.957    (.0253) 

o  = 

.05 

1.651   (.0509) 

1.673   (.0500) 

1.645   (.0515) 

a  = 

.10 

1.291   (.1013) 

1.312   (.1004) 

1.286   (.0993) 

.01 

112=13 

n2=14 

n2=15 

o  = 

2.306   (.0100) 

2.306   (.0101) 

2.318   (.0100) 

a  = 

.025 

1.954   (.0256) 

1.961    (.0249) 

1.958   (.0255) 

a  = 

.05 

1.642   (.0501) 

1.657   (.0501) 

1.665   (.0504) 

a  = 

.10 

1.297   (.1001) 

1.296    (.1005) 

1.305   (.1001) 
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